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Abstract 

A system modeling the electrophoretic motion of a charged rigid macro- 
molecule immersed in a incompressible ionized fluid is considered. The 
ionic concentration is governing by the Nernst-Planck equation coupled 
with the Poisson equation for the electrostatic potential, Navier-Stokes 
and Newtonian equations for the fluid and the macromolecule dynamics, 
respectively. A local in time existence result for suitable weak solutions 
is established, following the approach of [15] . 

1 Introduction 

Electrophoresis is the motion of charged colloidal particles or molecules through 
a solution under the action of an applied electrical field. This phenomenon is 
one of the most powerful analytical tools in colloidal science, being often used in 
the characterization of colloidal systems [33] . Modern applications include drug 
delivery and screening, manipulation of particles in micro-/nanofluidic systems, 
sequencing of genome of the organisms, forensic analysis, micro-chip design and 
others [8] , [32] , [20] . Despite the fact that the analytical description of this type 
of electrokinetic phenomenon is a difficulty task (due to the complicated inter- 
play between hydrodynamic and electric effects |41|). there is a vast literature 
on the subject. The majority dealing with special problems like electrophore- 
sis of bio-molecules with some geometric symmetry property or small surface 
potentials (see, for example, p], [19], [24], [35], [38]). Few of these papers deal 
with the rigorous mathematical aspects of the models involved. 
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In this paper we present a model of electrophoresis of a single particle im- 
mersed in a viscous and incompressible fluid (a ionic water solution). Unlike 
the models where the concentrations of charged particles range from colloidal to 
nano size (see [36]), we deal with the case where the colloidal particle is a very 
large rigid molecule and obeys the laws of classical mechanics. Also we consider 
the situation where the particle and fluid are contained in a bounded domain 
Oct 3 which represents the enclosure of the system (a tube, in the capillary 
electrophoresis). These are fairly reasonable assumptions in a number of prac- 
tical situations such us in the study of electrophoresis of proteins [T2] ■ On the 
other hand, we assume that the particle is far from the boundary of O. As a 
consequence we are able to analyze only the local motion of the particle, as we 
are considering the standard electrokinetics equations. In fact, in the standard 
models (see [41], p], [21]), the analysis is restrict to a single particle that is 
suspended in a fluid which fills all the region in M 3 exterior to the particle. This 
means that the boundaries are sufficiently far removed from the particle so that 
they have negligible effects on the electric and fluid velocity fields associated 
with the particle [23]. According to [37, when the geometry forces the mobile 
particle into close of proximity with a surface (which may itself be charged) 
complex phenomena arise and a more detailed analysis should be done. 

Although we are interesting mostly in the local motion of the particle, we do 
not neglect the inertial terms in the Navier-Stokes equations for fluid velocity 
and pressure. The Reynolds number for fluid flows in a typical electrophoretic 
motion is small and many authors consider the Stokes equations instead the 
Navier-Stokes equations (see [41], [32], [3])- However, it is not clear to us under 
what circumstances the linear theory can be used. With respect to this issue 
Allison and Stigter (2] suggest that the use of linearization of the constitutive 
equations with respect to the perturbing electric/flow fields are allowed provided 
the perturbing fields are weak. Also we believe that the analysis of the non 
simplified model considered here can be given us a support from the treatment 
of more general problems where the inertial effects could not be negligible (see 

m, 0). 

Another important aspect in the modeling of electrophoretic motion is the 
electrostatic potential. In more simplified models the total electrostatic poten- 
tial of the system is given by the applied external electric field (see [3], [38] ) 
and the local potential (due to particle charge distribution and ions) is neglect. 
On the other hand, the local electrostatic potential is usually modeled by the 
Poisson-Boltzman equation [5], [I], [7], [42]. This equation is derived from the 
assumption of thermodynamics equilibrium where the ionic distributions are 
not affected by fluid flows. According to [30], this is a reasonable assumption 
for the case of steady-state processes with stationary values of diffusive fluxes, 
but there are some important cases where convective transport of ions has non- 
trivial effects [27]. In this context, in the lines of [23] and [41], we consider 
a more convenient approach: the use of a convective-diffusion type equation, 
more precisely, the well known Nernst-Planck equation, for the ionic concentra- 
tion and Poisson equation for the total electric potential. Moreover, we only 
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impose C 2+a regularity on the domain occupied by the particle (without other 
geometric symmetry restrictions). 

Some of the mathematical aspects of a electrophoresis model based on the 
Poisson-Boltzmann theory have been discussed by the authors in the papers [4] 
and [5]. However, at least to the knowledge of the authors, there are no rigor- 
ous mathematical results on the Poisson-Nernst-Planck model of electrophoresis 
of the rigid macromolecules. From the mathematical point of view, the main 
feature relies in the highly coupled equations which must be solved in time de- 
pendent domains. We establish the proof of a existence theorem (see Theorem 
13.11 in Section 3) for appropriate weak solutions for this system. The proof is 
based on the approximation technique introduced in the references |14j , |15j for 
the study of the motion of rigid bodies in viscous fluid. A sequence of approxi- 
mate smooth solutions is construct and the solution is obtained as the limit of 
this sequence. Our scheme is local in time in the sense that the macromolecule 
does not touch the boundary of the enclosure in the time interval of existence. 
Employing a linearization and regularization procedure, this permits us to ob- 
tain suitable energy bounds for the electric potential, fluid velocity and ionic 
concentration (see Lemmas 14.21 14.31 14.5[) . As a consequence, we establish spe- 
cial results on time compactness (see Lemmas 15.11 and I5.2[) that give us strong 
convergence results for fluid velocity and ionic concentration. Then, the main 
result follows from the passage to the limit in the approximation equations (see 
Lemmas 15. 3[ 15.41 and 15. 5|) and is announced in Theorem 1 . 

The outline of the paper is as follows. In Section 2 we introduce the cou- 
pled system of governing equations; in the Section 3 we obtain formal energy 
estimates that give us a weak formulation for the system; in the Section 4 we 
construct a sequence of approximate solutions that are uniformly bounded in 
energy spaces; in the Section 5 we obtain appropriate compactness results that 
permit us to pass the limit the approximate solutions. 

In a subsequent work we will investigate the effect of proximal boundaries in 
the electrophoresis of many particles and consider the possible contact between 
two particles or with the boundary, in the lines of |18j . 

2 Governing Equations 

We assume that the solvent and macromolecule occupy a bounded open con- 
nected set 0cK 3 and the solvent contains J ionic species. At the initial time, 
the macromolecule occupies a compact region Kq, where Kq CC O is a open 
connected set such that 

70 :=dist(K ,dO) > 0. (2.1) 

The fluid domain at the initial time is denoted by fio := 0\Kq. We also assume 
that Kq and fta are C 2+a -domains, < a < 1, and define K t = Q(t)Ko as the 
position of the particle at time t, Q{t) is an affine isometry such that Q(0) = /. 
Also, let us denote = 0\K t , i.e., the fluid domain at time t and assume 
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(formally) that there exists T > such that 

7 (t) := fist(K u 00) > 7, < 7 < 70 (2.2) 

for all t S [0,T]. 

It is important to remark that, although we have not imposed any sym- 
metry hypothesis on Kg, the regularity imposed is not satisfactory in a large 
number of situations. In fact, the determination of the electrostatic interface in 
biological processes is not a trivial task. According to pj] (see also [12]), any 
model of molecular surface must follow, in some way, the boundary of the space 
from which other molecules are excluded. Usually, the molecular surface (the 
accessibility surface of the solvent) is taken to be the surface described by a 
point on the surface of an idealized spherical solvent probe as it is rolled around 
the solvated molecule. This lead us to a C° surface with the eventual presence 
of cusps and with multiple connected components. Accordingly, we are dealing 
with a idealized situation in order to obtain appropriate regularity results and 
a priori estimates for the approximate solutions of the system. However, we be- 
lieve that this condition can be relaxed using a smoothness domain technique. 
We will discuss this issue in a subsequent paper. 

If the dielectric constant of the particle and of the solvent are K\ > and 
k 2 > 0, respectively, the electrostatic potential ip, for each t £ [0, T], is governed 
by the Poisson equation 

,/ 

V • (k(x,t)VV>(x,i)) = -Aire V^(x,t) - 4np(x.,t), x e O, (2.3) 



where 



*(x,i) = { ^f^*^' (2.4) 
0, x € K t ; 



k:Ox[0,T]-> L(R 3 ,R 3 ) is given by %(x,t) = <%/s(x, t), with 

«(x,t) := Sij{KiX Kt + K 2 Io t )(x,<); 

e is the electron charge, Zi and Mi are the valence and concentration of the zth 
ionic species, respectively. Moreover, 

P ( x ,t) = { MQ-H^xe^, 

\ o, xefi*, y ' 

where po is a fixed charge distribution on K such that supppo CC K . Note 
that (|2.5|) implies that the fixed charges of the particle are invariant under rigid 
motion. As a consequence, we have the conservation of total fixed charges 



/ pdx ) (*) = / p (x)dx, 
Jo / Jo 



(2.6) 



for all t G [0, T]. These assumptions on p corresponds, of course, to an idealized 
situation. In fact, the determination of the properties of fixed charges in colloidal 
molecules is a very difficulty task (see [12], |28|). 
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We assume transmission boundary conditions for ip 

M^t) =V»i(x,f), xeM,, (2.7) 

{Kidu^i - K 2 d^ 2 )(x, t) = 0, x e dK t , (2.8) 

$j(x,i) = *(x), xe30, (2.9) 



where is the exterior normal to dK t and ^ represents a stationary elec- 

trostatic potential on dO, ipi = ^|x t j V*2 = ''/'Inf 

The evolution of A/i, for each z G {1, . . . , J}, is described by the Nernst- 
Planck equation 

dtMfrt) = V • (-NiVf + dkVNi + ^jVivJ (x,i), (x,t) G fit x (0,T), 

(2.10) 

where dj is the ionic diffusion coefficient of the ith ionic species, Vf is the 
solvent velocity field, Kb is the Boltzmann constant and 9 is the temperature 
of the system (which we suppose constant). It is important to remark that, 
for obvious physical reasons, we seek non-negative functions Mi in fi t x [0,T]. 
The boundary conditions correspond to no ion diffusion and no ion conduction 
through the boundaries (see [23]) are given by 

(duAfi + —M % dyih) (x,t) - 0, (x,f) G (dO U dK t ) x (0,T). (2.11) 

The system is complemented with the initial condition 

A/";(x,0) = M,o(x) > 0, a. e. x G fi . (2.12) 

The solvent velocity field is governed by the Navier-Stokes equations for incom- 
pressible flows 

M/(ftv / + V-(v / ®v / )(x,i)+ 

-r ? Av / (x,t)+VKx,t) =/I/F(x,t), (x, t) e fit x (0, T) (2.13) 
V-v/(x,t) =0, (x,t) G fit x (0,T), 

where 

F(x,t) = -5^e(g i V^ 2 )(x,t) (2.14) 

i=l 

is the electrical forcing term, p is the pressure, r\ > and /I/ > are the viscosity 
and the mass density of the fluid, respectively. 

Let denote by x c (t), v c (t) and w(t) the center of mass, the translational 
velocity and the rotation vector of the particle in the time t, respectively. If A, 
v p and ]I p > are the symmetric inertial matrix, velocity and the mass density 
of the particle, we have 

y T Ay=JI p f |y x (x-x c (0))| 2 dx, Vy G K 3 , (2.15) 

and 



v p (x,i) =v c (*)+w(i) X (x-x c (i)), (x, t) £ K t x (0, T). (2.16) 
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Using (j2.3[) and (j2.14j) . the electrical forcing term F can be written in the di- 
vergence form P = V • <je, where <te — ((&E)ij) is the well known electrostatic 
tensor, 

(see [H]). If M is the mass of the particle, the evolution law for its motion is 
given by the Newtonian dynamic equations 

M^-(t)=[ <7fr(x,t) ■ f(x,t)ds(x) + [ a E {x,t) ■ is{xi,t)ds(x) (2.17) 



and 

A 

dt 



A-^-(t) = I (x-x c (t)) x (a H ■ v)(x,t)ds{x) +w x (A-w)+ 

dKt (2.18) 
(x - x c (t)) x (cr B • z/)(x,t)ds(x), 

where <jh is the stress tensor of the fluid. If we set D(vf) — |(Vv/ + (Vv/) r ), 
then 

a ff (x,t) = 2rjD(v/(x,t)) -p(x,t)L 

We assume the homogeneous Dirichlet condition for vj on dO and require 
the velocity and normal stress to be continuous at the interface between the 
solid body and fluid 

v/(x,t) = 0, (x,t) e 90 x (0,T), (2.19) 
v/(x,i) = v p (x,t)eM f x(0,T), (2.20) 
(K + a £ )-F)(x,() = (E^)(x,i),(x,()eSJ( t x(0,r), (2.21) 

where E is the Cauchy stress tensor in the solid. Also, we have the following 
initial conditions for the velocities 

v c (0) = v c , , w(0) = w , 
v p (x, 0) = v p ,o(x) := v c , + w x (x - x c (0)), x £ K , ^ 
v /( x , 0) = v /,o(x), x6!J , 
v/ )0 (x) = v p ,o(x), x e dK . 

In the "zeta" potential approach for electrophoresis modeling, (|2 .20[) is re- 
placed by a "slip" boundary condition on vj: a nonlinear Dirichlet condition 
that depends on ip and it is based on a Prandtl boundary layer approximation 
(see [3], [4]). This approximation is relatively accurate when the Debye screen- 
ing length of the macromolecule is much smaller comparable with its radius of 
curvature |38) . Otherwise, in the more realistic cases, this approach is no longer 
correct l23l. 



3 Weak Formulation and Main Result 



We can obtain a weak formulation for the problem (|2.3[) - (|2.22j) if we take into 
account the energy framework of the system. First, we need to introduce a global 
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formulation for (|2~T3)) . f2~T7|) - (j2~T^)) . Following [H], we define in O x [0,T] the 
Eulerian densities 

H P (x,t) = ~p, p T Kt (x,t), = 7*/^n t (x,*), 

the global density /i = /x p + and the global velocity 



u(x,*) = 



v p (x, t) if xeK t , 
v/(x, t) if x £ (V 



In view of mass conservation, /i must satisfies 

9tAt + V • (jm) = 0, (3.23) 

in V'(0 x (0,T)). Also, in terms of the Eulerian quantities, (|2~T3)l . (|2~T71) and 
(j2~18l) can be expressed, in V'{0 x (0,T)) 3 , as 

dt(iifu) + V • (/i/u ® u) = ^-V • (fJ>fcr H ) + rr£ • V^ p + V • (^fa E ), 

M/ Mp 

dt(Lipu) + V • (/x p u® u) = —V • (/ipS) - — cr/f • V/x/ - — (tb ■ V/i/, 

M P M/ M/ 

where we have consider the balance the momentum in the fluid and particle 
(see [S]). Summing the above equations and using (|2.2ip we obtain the global 
formulation for (j2~T3l) . (|2~T7)) - ([2~T8)) . 

dt(fm) + V • (/iu ® u) = V • (^rHfa H + + M/CTb J (3.24) 

V A 1 / M P / 

in P'(0 x (0, T)) 3 . Furthermore, from ([2~2Uj) . (|2"36]l , f2~T^)) and (|2~2"2"|) it is clear 
that, in the sense of distributions in O x (0,T), 

Vu = 0, /vD(u)=0, (3.25) 

and that 

u(x, 0) = u (x) := (l Ko v pfi ){x) + (Xo v /i0 )(x), x e O 

/x(x,0) = ^o(x) :=-p f X Ko (x) +7J p 2n (x), x e O, (3.26) 
u = in aOx(0,T). 

Formally, taking the inner product of (|3.24j) with u, integrating by parts, 
using (|3l>5]) and (j3~23| we obtain, for all t € (0, T), 



(3.27) 



(^J n f Du : Dudx^j (t) - fi f F ■ udx^j (i) = 0. 



Note that E = E k + E d + E p , where 
E k \p,,u,Ni,...,Nj,i>]{ t ) 



\ / o mM 2 ^) (*) 



= Q^ M/ |u| 2 dx) (t) + i(M|u c (t)| 2 +w(if„4w(t)) 
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is the kinetic energy at time t, 

_ r 

is the viscous dissipation and 



£ rf [/i,u,A/i,...,A/j,V>](£) = [ [ HfDu : Dudxdr 

M/ Jo jo 



E p [fi,u,N'i ) ...,J^j,ip](t) = - / jti/P-udx 

Jo Jo 

is the total work of the electrical force. Then (|3. 271) corresponds to the energy 
conservation of the system. Moreover, integrating (|3.27|) with respect to t we 
obtain the energy (mechanical) bound 

1 / /i|u| 2 dx) (t) + r) [ [ |Vu| 2 dx < f ( (ifF ■ udxdr + i / ^ |u | 2 dx. 

2 Jo / Jo Jo Jo Jo ^ J 

(3.28) 

As a consequence, recalling that there exists a positive constant C* = C*(0) 
such that 

||u(.,i)|[ 0)2 ,O < C*||Vu(.,i)[[ , 2 ,o, 

it is easy to obtain the estimate 



/ |u| 2 dx)(i)+/ / |Vu| 2 dx<%^/ ||F(.,r)|| 2 , 2inT dr 

JO / JO JO Umin JO 



(3.29) 



|llo| MX, 

o 



where fi min = xain{2r], /i p , fi f } and j/ moa! = max{/i p , 
Now, supposing that * G H 1 (dO), we define 

= {$ g ff^O), $|ao = *}. (3.30) 

If we take the product of (|2.3[) with ip — iff, where g L*, we obtain, using 
(I2~7l) and (1231) . 



y k|VV>| 2 cZx^ (t) -4neJ2 {^J ZiKi>dxJ (t) - 4tt (^J pijjdx^j (t) 
= -4tt (^J p^dx^j (t) - 4:TreJ2 {^J ZiNi^dxj (t) + nVij} ■ V$dx^j (t), 
for each t G [0, T]. Using Cauchy and Young inequalities, we have 
i ^ K |W>| 2 dx^ (i) -47re^ ^ Z^d^j (t) + 
-4tt ^ pi/;dx^ (t) < i K |V$| 2 dx^ (3.31) 

~ 4?re (/ z ^* dx ) (*) ~ 4?r (/ /°* dx ) (*). 
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for each t € [0, T], The left side of (|3.3ip corresponds to the electrostatic energy 
contribution (the negative of the free energy according to [33]) which we denote 
by E e i\ the first term in E e i corresponds to the self-energy of the electric field 
and the next two terms are the electrostatic energies of the ions and fixed charges 
(see [T2] , [25] , [33] ) . As a consequence of (|3.31[) , for known A/i and u and each 
t 6 [0,T], 

E rf [M,u,M,...,JV},V](t) < S e; [/i,u,M,...,A0,$]W, (3.32) 
for all §eL t . 

Note that (|2.10[) (as well as (13.231) ) corresponds to a conservation type equa- 
tion. In fact, from the transport theorem 

j t Qf Nidxj (t) = Qf 9 t A/idx^ (t) + Qf^ Mv/ • i/ds) (t), 

then, integrating (|2.10[) in O t x [0,t], using Gauss theorem, (|2.1ip . (|2.19[) and 
(|2.12j) . wc obtain the molar conservation for the ith ionic species, 

( f Mdx) (i) = / Nifidx, Vt G [0,T]. (3.33) 

In particular, using (|2.6|l . we have the total charge conservation for the system. 

From the above discussion, in particular, from (|3.32p and (I3.27[) , we have 
a natural weak formulation of the above system (see Definition 13.11 below), 
which is obtained from the following minimization problem with constraints: 
Find (/i, u, A/i, . . . , A/j, tp) ( m appropriate functional spaces) that minimizes the 
energy E in O x (0, T) and such that minimizes E e i [//, u, A/i, . . . , Afj, .](t) 

in L*, for each t G [0,T]; (j2~T0t . (|2~TT]) - (|2~T2|) . (j3~25]) - (j3~26| are the constraints; 
furthermore, the pressure p and the Cauchy stress tensor E are the Lagrange 
multipliers of the problem. 

The additional hypotheses on the data are 

H po e L 2 (0), supppo CC X , * G ff 1 ^©), v /l0 G L 2 (n ) 3 , v L0 \ao = 
M,o G L 2 (0), suppM,o CCflo, A/"i,o > a. e. in Q . 

Despite the fact that the ions are concentrated very close to dK t , the hypotheses 
on the support of po and A/^o are reasonable as there is an ion exclusion region 
close to the boundary of the molecule [12] . 

In order to obtain the suitable functional spaces of the solutions, we observe 
that (|3.28|) and (|3.29[) can give us if ^bounds for u and ifj- However, we need 
uniform L 2 -estimates for F and for Mi- As we will see, for small T (depending 
on the initial data) uniform L (r2 t )-estimates for Vi/> and A/i can be obtained; 
as a consequence, we can estimate F and the conduction term in the weak 
formulation for (|2.10[) . 

We need to define the following functional spaces 

L\Q,T;H\nt)) = L 2 (h T ), Vv £ L 2 (f7 T ) 3 }, 

L°°(Q,T;LP{n t )) ={vG L p (^t), ess sup \\v(.,t)\\ . p .n t < oo}, 

te(o,T) 
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where 

n T = |J {<} x fl t . 

te(o,T) 

Definition 3.1. Let us fo assume that Kq and f2o are C 2+a -domains, < a < 
1, p is invariant under rigid motion (see h2.5\) ). hypothesis H and V2. 1)) . Then 
(p., u,^,M, • • • , AO) is a weafc solution of K3.23]) . {3.21$ , \3.25)) - \3.2G\) . coupled 
with PTTP- iTOj) . fOfl) . WJT\) - WJm in (0,T) i/if safe/es 

L (leL^Ox (0, T)), u G L°°(0, T; i 2 (C)) 3 n L 2 (0, T; ^(C)) 3 , 
V- G L°°(0,T;ff !((!?)), w> g L°°(o,r;L 4 (r! t )) 3 , 
M G i 2 (0, T; F x (n t )) n L°°(0, T; L 2 (tt t )), V* G {1, . . . , J}. 

2. For F given in \2.11$ , (u, fj,) satisfies 

n(/xu ■ 9jW + /iu ® u : D(w) — T]D(u) : D(w) + /z/F • w) dxdr 
j 

+ J /i Uo • w dx = ^ ^u-wdx^j (t), 

(3.34) 

for all w G 5 and a.e. i G (0,T), where 

S = {w G fl^O x (0,T)) 3 /w(.,t) G S(t), Vt G (0,T)} 
and S(i) = {w G iJ 1 (O) 3 /V • w = 0, fi p D(w) = 0}. Furthermore, 

dt/j, + V • (/xu) =0, V • u = in D'(Ox(0, T)), (3.35) 
/x p £>(u) = in V'(0 x (0,T)) 3 . (3.36) 

3. For a. e. t £ (0,T), ip(.,t) is the solution of the 

(^J K|VV>| 2 dx^ (t) -47re^ (^j Z.M^dx^j (t)+ 

- 4tt (^J pi>dxj (t) = -Air (^J p^dx^j (t)+ (3.37) 

- 47re^ ^ ZiAT^dxj (t) + (^J K^tp ■ Wdx^j (i), 

for all §6Lt={f e H l (0), $| ao = 
4- For each i G {1, . . . , J}, Afi satisfies 

( [ Afiqdx) (t) + [ [ ? (u • VM)dxdr+ 
\Jn f / Jo Jn T 

+ di L I ( VA/ " 4 + ?e UlV ^ ' Vq dxdT+ ^ 3 ' 38 ^ 

nAfi q dxdr — I / A/i,o To^x ) 

/or a.e. i G (0,T) and /or aZ/ ^ G if (fiy). Furthermore, Afi > a.e. in 
fi T and (^f nt Af 2 dxj (t) = J^Af^odx, for a.e. t G (0,T). 
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The main result of the paper is given below 

Theorem 3.1. There exists T* > and a weak solution of h3.2']\) - ^3~Wjl , cou- 
pled with $EW-(MIM), W^- W^ in (0,T*). Moreover, 

T* = ia£{t > 0;dM(K t ,dO) > 0}. 

As remarked in the introduction of the present work, the strategy of the 
proof follows the lines of the results in references [14] and [15] for the mo- 
tion of rigid bodies in incompressible fluids. It is based on the solution of 
an approximated system, which give us an approximate sequence of solutions 
^„,A/"i,„,...,A0,n) for (13341) - (13381) . The solution (u, n, tp, Aft, . . . , Nj) 
is built up as a limit of these approximations. 

The new feature here is the coupling between (|3.34|l - (|3.38|l which lead us to 
consider an appropriate fixed-point schema in order to construct the approxi- 
mate sequence. The existence results for the Lagrangian version of the parabolic 
problem (|3.38|) depend on the suitable estimates on the term containing VV>2 
(see Lemma I4.3[) . Then, the crucial point of the proof is the obtention of the 
uniform L 4 -estimates for Vip2 (see Lemma [4.11 below), in the sense that the 
bound does not depend on the particle motion. Using results on singular inte- 
gral operator and elliptic regularity, this is possible if we assume that (]2.2[) is 
valid. 

4 Approximate solutions 

Let us fix < 7 < 7o- We begin with a technical result 

Lemma 4.1. Suppose t > and Q(t) an affine isometry such that Q(0) = I 
and K t — Q(t)Ko, Q t = 0\K t satisfy V2.2]) . Consider the situation in which 

J 

-Ane ^2 1i ( x ; *) - 47rp(x, t) 
i=i 

is replaced by f(-,t) G L 2 (0) and ^ is replaced by g G H 1 (dO). Then the 
problem \3.3 r ty has a unique solution ip(.,t) G i? 1 (C) that satisfies 

max{||VV'i(.,i)llo ) 4,*«> ||VV> 2 (.,i)lkW < C l (||/(.,t)||o,a,o + ||ff||i,2,8o), 

(4.39) 

where d = C 1 {O,K ,Ki,K 2 ,i). If, f(.,t) G C a (0) and g G C 2+a (dO), we 
have 

(M-,t),M,t)) G C 2+a (K t ) x C 2+a (Q~t). (4.40) 

Proof. From a theorem due to Stampacchia (see [6]) it is a routine to check 
that there exists an unique solution in iJ 1 (C) for the problem (|3.37p . Let us to 
extend / to be zero outside of O. Then, if w(x, t) = —(G * /)(x, t) where G is 
the fundamental solution of the Laplace equation in M 3 , G H 2 (0), 

Aw(x,f) = /(x,t), a. e. x 6 M 3 , (4.41) 
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and there exists C = C(O) such that 

IM.,i)lk2,o<C||/(.,i)||o,2,o. (4.42) 

For the above results see Theorem 9.9 in [21] . 

We define B 7 = {xeO,0< dist(x, dO) < 7/2} and consider the following 
problem 

A$(x,t)=0, x£S 7 , 
$(x,i)=0, x€dB^\dO, (4.43) 
$(x, i) = K 2 g(x) - k 2 " 1 oj(x, t), x 6 90. 

From well known results (see Theorem 2.2 in Chapter 4 in [29 ), this problem 
has a unique solution 

$(.,t) e H 3/2 (B 7 ) (4.44) 

that satisfies 

ll*(.,*)ll3/2 A o < C(\\g\U,2,oo + M.,t)\\i,2,do), (4.45) 

where C — C(0, K2, 7) and we have extended $ to be zero in 0\B 1 . 
Now, let us to consider < a < 7/4 and define the set 

A a {t) = {x e fi t , < dist(x, aXt) < cr}. 

Note that, from (|2.2[) . we have dist(Ar(*), B 7 ) > 7/4. We also define 

and consider the following auxiliary problem 

A0i(£,i)=O, £etf , 

( K 2 1 ^ 2 -Kr 1 0i)(^t) = (Kr 1 -«2 1 M^*). ^ea^o, (4.46) 
(9^ 2 -^ 1 )(e,t) = o, ee^o, 

&(£,*) = <), ee9Ar(0)Wo. 

This problem has a unique solution 

(M,t),M,t)) e ff 3/2 (^o) x ff 3 / 2 (E 3 \Xo), (4.47) 



where we have extended to be zero outside A a (0) U JsTo- In fact this follows 
from suitable modifications in the arguments of the Theorem 7.2 in the reference 
[40] (see also [22]). Furthermore, from the results on singular integral operators 
of [TU] and [10] (for related results see [1]), there exists a constant C > 0, 
depending on Ko, K\, K2 and 7 such that 

\\M;t)h/2.2,K < C||S(.,t) || 1,2,9^0, 

(4.48) 

ll^2(.,*)|| 3 /2,2,K3\K'o - ^II^Gj^lll.a.SKo- 
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Now, if we define </>;(x,i) = 4>i(Q 1 (i)x, t), i = 1,2, we see that 

A&(x,t)=0, xG-fQ, 
A</. 2 (x,i) =0, xeA CT (t), 
(«^V2 - «rVl)(x, i) = (k^ 1 - /s^ 1 )w(x 1 i), x G 0Jf t , (4.49) 
(9^ 2 -5^!)(x,t)=0, xe9if f , 

02(x,t) = O, xe9i ff (f)\M f . 

This follows from an elementary calculation using the properties of Q(t). We 
claim that 

ib(xt) = i V'i(x,i)=% 1 (^+^ + *)(x,*), xGJT t , 
\ V 2 (x,t) = ^ 1 (0 2 +w + $)(x,t), XGO t . 

In fact, the transmission boundary conditions (|2.7|) - (|2.9p are satisfied in the 
sense of traces, in particular, (|2.7p implies that i/j(.,t) G also it it easy 

to check that (|3.37[) is valid. The result then follows from the uniqueness of the 
problem (|3.37p . Now, note that, 

VfS5(e, t) • V c w(C, t) = V x w(x, i) • V y w(y, t), (4.51) 

for x = Q(t)£, y = Q(t)C> V£, C G 2£o- Then, if we consider the Slobodetskii 
norm 



l|2(.,t)||l/ a , a ,Xo = / 2 2 (C,t)de 



/ |v c 2(c,t)i 2 de+ 



Kb 

if •'Ko Is — Si 

it is easy to check that 

l|w(.,t)||! /aj a,* = I|w(.,t)ll3/2,2,jc t . ( 4 -52) 

Henceforth, from (|4.5ip . standard Sobolev embedding, (|4.48l) . the trace theorem, 
(|432|) and (|4~42| . we have 

(f _\vM^t)\ 4 dt) =([ _|v € &(£,f)|W) 

\jR 3 \K t ) \JW\K a ) 

< C\\M.,t)\\ 3/22 ^ x -^ < C||w(.,t)lka,e* ( 4 - 53 ) 
<C , ||w(. ) t)||3/a,a,jr =C f |K J t)||3/2,a,ic t 

<C||w(. ) t)||3/ 3 ,a,O<C||/(.,t)|| ,2,O, 

where C = C(K , O). Similarly, 

\\VM-Ml A ,K t <C\\f{.M^o- (4.54) 

Consequently, (|4.39[) follows from (|4.50p . (|4.45p . the trace theorem, Sobolev 
embedding, (|4~42]) . (|4"33j) and (|4~54| . 
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Now, let us to suppose that f(.,t) G C a (0) and consider the extension of 
/(.,*) to M 3 (until denoted by /) such that f(.,t) G C£(M 3 ). Then, we have 
uj(.,t) G C 2+a {0) (see Lemma 4.4 in [21]). If g G C 2+a (dO), from classical 
results, 

(M-,t),M-,t)) g c 2+Q (7^) X c 2+a (M 3 \K t ) 

(see Theorem 2.1, Chapter 14 in [26 ) and 

*(.,t) g c 2+a (K;) 

(see Theorem 6.14 in [21]). Hence (|4.40|) follows from (|4.50|) . □ 

In what follows we describe how to construct the sequence of approximations. 
We adapt the proof in [15]. The idea is to introduce an approximation scheme 
which consists in solving a system of regularized equations: an "almost" linear 
problem related to the velocity field as well as appropriate linear problems for 
Afi and ip. The existence of these regularized solutions is obtained from the use 
of a fixed point type theorem. As described below, this is done in small time 
intervals, chosen in a such way that the advecting vector field is close to the 
identity, (12.2[) is valid in each time interval, and using a Lagrangian Galerkin 
type method for the linear problems related to u and Afi (a similar approach 
may be seen in the study of the free surface problems [39]). 

Let us consider 

L 2 > ~T7T max {||M.,o||o,2,n , • • • , ||Ajo||o,2,n }) ( 4 - 55 ) 
d 1 

mtn 

where d m i n = min{l/2, di/8, . . . , dj/8} and 

A = Z2e 2 d Z 2 max C 2 ^ + llp4h ^ + ml2 dol (4 56) 

where C2 = 2C 2 max{327r 2 e 2 J^/ =1 |Zj| 2 , 32-7T 2 }, d max = max{di, . . . , dj} and 
Z m ax = max{|Zi|, . . . , \Zj\}. 

From the Sobolev embedding theorem, there exists a positive constant C„» = 
C**(r2o) such that 

H/llcMA, <a*||/||i,2,cv (4-57) 
Let us define L\ > such that 

Ll > ^clc^LlTi^Zi? {L% + + ml2 do)+ 

+ f | Uo |2 dXj 

Mmin JO 

and 

0<T<min{^,JL}. (4.59) 

Following [T5] we observe that it is possible to show the existence of a home- 
omorphism from the space of incompressible vector fields 

u G L 2 {0, T; H^{0)) 3 n L°°(0, T; L 2 (0)) 3 



15 



and which corresponds to a rigid motions in K t into the representation space 

Y T = C°' 1 ([0,T];M 3 xR 3 )x{v £ L°°(0, T; L 2 (fl a )) 3 nL 2 (0, T; Hq(£1 )) 3 , V-v = 0}. 

Let us also define Wt = L°°(0, T; L 2 (flo)). Here we will consider the natu- 
ral corresponding norms in Yt and Wt- Now, using standard regularization 
techniques, we consider 

• Uo. e G C^°(0) 3 such that V • Uo, c = in O and corresponds to a rigid 
motion in Kq\ furthermore, Uo, e converges to Uo in L 2 (0) 3 as e ^ + . 

• M,0,e 6 C 2+a (0), suppM,o,£ CCllo, M,o, e > in fi , M,o,e converges to 
M,o in L 2 (C) as e ->■ 0+. 

• * £ G C 2+Q (9C), * e converges to * in H\dO) as e -» 0+. 

• po,e € Co (-^o)j Po,e converges to p in L 2 (K ) as e — > + . 

Also we consider, for each / G L°°(0, T; L 2 (0)), r e (/) G L°°(0, T; C Q (0)) such 
that 

|ke(/)|U~(0,T;i2(O)) < ||/||z,°°(O,T;Z, 2 (0)), ( 4 - 60 ) 

for e > sufficiently small and r e (/ e ) converges to / in L 2 (0 x (0,T)) with / £ 
converging to / in L 2 (0 x (0, T)) as e — > + . 
Let us also consider, for / G L 2 (0,T), 

[fY(t)= / T g e (t-r)f(r)dr, (4.61) 



where {<? e } c >o is a family of regularizing kernels with respect to time such that 
[/] e (.) G C°°[0,T]. Furthermore, if / G i 4 (0,T) 

II [fY \WW) < ll/llo,4,(o,T) (4.62) 

for e > sufficiently small and [f e ] e converges to [/] in L 2 (0,T) with f e con- 
verging to / in L 2 (0,T). 
Now, as in [15], for any 

v = (x c (i),0(*),v)GF T , 

we consider the incompressible field 

v = 6(v) in L°°(0, T; L 2 {0)) 3 n i 2 (0, T; ^(O)) 3 

and its regularized version v e = JZ € (v), where lZ e is a regularization opera- 
tor (see [H]), such that TZ e {v) is analytical in time and smooth in space (in 
particular, Lipschitz in space), lZ e (v e ) — >• v in i 2 (0 x (0,T)) 3 if v e — > v in 
L 2 (C x (0,T)) 3 , as e ^ + . Furthermore, v e is divergence free in O and corre- 
sponds to a rigid motion in the particle domains K Ci t — {X e (£,0,t), £ G i^o}, 
where X e is the Lagrangian flow of v e : for each < s < T and £ G 0, 

4y e (£,s,t)=v e (X e (£,8,t),t), 0<t<T, t^s 

(4.63) 
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Also, we denote, for each t € [0,T], Q, e j — 0\K e ^ t - From the construction of 
(see [15]), if v = 0(v) for some v £ Yt and ||v||y T < L\, we have 

ll v e||L^(0,T;L 2 (C>)) 3 nL 2 (0,T;Hi(C>)) 3 < ^1 (4.64) 

for e > sufficiently small. Then, let consider e' > such that (|4.60j) . (|4.64p . 
(|4^]) and 

||Uo, e ||o,2,0 < ||u ||o,2,e>, ||M,0,e||o,2,0 < ||M,o||o,2,0 

(4.65) 

||*e||l,2,flO < 11*111,2,90, ||po,e||o,2,0 < ||Po||o,2,0 

are valid for all < e < e' . As a consequence of (|4.64[) and (|4.59[) . it is easy to 
check that 

7e (f) := dist(dK e , t , dO) > 7 (4.66) 

for all t £ [0,T] and < e < e'. 

Now, for < e < e' fixed, let us define F e = VA e . From the smoothness of 
v e , (|4.64p and (|4.63p . it is possible to choose N £ N* (depending only on e and 
Li) and t = T/N such that 

sup sup ||I-F e (.,mto,t)|U- (0) < 1/3. (4.67) 

0<m<N t£[mt ,(m+l)t ] 

The inequality (|4.67[) implies that, in each time interval (ifc-i, tk), the advecting 
vector field is close to the identity. This is an important point in the proof 
and it is based on [!!>]■ As we will see, from (|4.67[) . the Lagrangian forms 
of certain operators that appear in the Lagrangian formulation for (|3.37[) and 
(|3.38|) are uniformly elliptic. In fact, from (I4.67|) . it is easy to check that, for 
each m £ {0, 1, ... , N} and t £ [mto, (m + l)io], 

l/4|y| 2 < iF^^mto.tJ-yl 2 < 4|y| 2 , a. e. (eO, (4.68) 

Vy £l 3 . As a consequence, we have existence results for the related linearized 
problems, in each time interval 

Let us consider L2, Li, T, e' > 0, N and to a s above. For each < e < e' 
and in each time interval (tk— lj tfe), we want to solve a set of linearized problems 
and to apply the Schauder's fixed point theorem. Below we give the details for 
the first time interval (0,io). 

We take (v, $1 , . . . , 1? j) £ B x X x ... x X, where 

8 = {ve Y to /\\v\\ Yto < LJ, X = {di £ W t J\\Uw t0 < L 2 }. 

Extend v(.,t) to be zero for t £ (to,T), we want to solve (in the given order) 
the following linearized problems 

PI Problem ([5371) in (0,t ) with Af { replaced by 



0,, e (.,t) :=r e (0i)(.,t), 



(4.69) 
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where = "3i(X e (.,t, 0),t) and di(.,t) is extended to be zero outside 

of fio! the domains (£2$, At) replaced by (0 £)t , -f£T e ,t), ^ replaced by \l/ e , po 
replaced by po,£i K replaced by 

K e (x, t) = n{lK t . t (x, i) + K2Zn e _ t (x, i). (4.70) 

and p, qi replaced by p e , Qi )£ give as in (|2.5|) and (|2.4[) (with the suitable 
modifications), respectively. Let us to denote this solution as ijj e . 

P2 Problem (j33gj) in (0,f ), with (-0,u) replaced by ([-0 e ] £ ,v e ), 7V ii0 by A/"i,o,e 
and the domains (fl t ,K t ) replaced by (fl e ,t, Kt.t)- Let us denote this 
solution as A/i, e . 

P3 Problem (pTM)) in (0,t ) with ^(x,i) replaced by 

Me(x,t) =/x 1 ljf 6t (x,t) +~p 2 Xn^ t (x ) t) (4.71) 

(the solution of (|3.23[) . considering v e instead u), F by F e , Uo by Uo, e and 
replacing the term /ju <g) u by ^ £ v e (g> u. Here F e corresponds to F given 
in (j2~Tl|) with (-0,M) replaced by (if e ,jV ¥ ). 

4.1 The Poisson System 

Lemma 4.2. T/ie problem PI /ias a unique solution 6 _L°°(0, io; -^(O)) 
</iaf satisfies 

(iM-> t), V>2,e(., t)) e C 2+q (1?m) x C 2+a (T!^), (4.72) 

/or a. e. i e (0,io)- Furthermore, for a. e. t £ (0,io) 

max{||VV 1 , e (. ) t)||^4,K et ,||VV2 J e(.^)llo,4 > n et } 

(4.73) 

<C 2 (L 2 + || /90 || 2 !2 ^ + ||vI/|| 2 !29o ), 
w/iere C 2 = 2C 2 max{327r 2 e 2 £f =1 | ^ | 2 , 327r 2 e 2 }. 

Proof. First we observe that, for all £ G [0, to], we have V • v e (.,t) = 0, so that 
det F e (., t, 0) — 1, as is easy to check. As a consequence, from (|4.69[) and (|4.60|) . 
we have 

||<M->*)llo,2,n e . t <L 2 (4.74) 

for a. e. t £ (0,to). Furthermore, from the smoothness of X e and the definition 
of r e (tii), tf;, £ (.,t) £ C a (n^), a. e. te (0,t ). Then, recalling fil^o) . the results 
follow from Lemma |4. 11 with 

./ 

/ = -47re^ ?i) e(.,t) -47rp e (.,t) 

i=l 

and 5 = >J/ e . In (|4.73|) we have also used (|4.65j) . □ 
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Remark 4.1. We have also the estimates, for a. e. t G (0,to), 

||cJ £ (.,t)|| C 2 + « ( o) < Ce(^2 + ||po||o,2,K ) 5 

Pe(.,*)|lc2+«.(ZJ) < C(||* e || C 2 +Q(90) + ||w e (.,t)||c. 3+ a (SO) ) 

< C € {L 2 + ||po||o,a,Jt + ll*l|i,2,ao), 

||^2,e(-,t)Hc=+*(a^) ^ C||^e(-^)||c2+<»(aiC M ) < C 6 (L 2 + ||po||o,2,Ko), 

where C e = C e (0, Ko,e, Ki, K2,J, Z±, . . . , Zj). As a consequence, for a. e. 
* 6(0, to), 

HlM.,*)llc» + « TO < C^ 1 + IIpoIIoaKo + ||*l|i,a,ao), (4.75) 

where C e = C £ (0, Ko, e, m, «2,7, £2, Z\, ■ ■ ■ , Zj). Moreover, for each £ G O 
the function ^ e (£, .) is measurable with respect to the variable t £ [0, to]. The 
same is valid for the first and second derivatives of ipe(£, ■) with respect to £. In 
fact, this is a consequence of (|4.50|) as w e (£, .) is clearly measurable and $ e (£, •), 
4>e(£> •) can be represented in terms of simple and double layer potentials with 
measurable time dependent density functions (see the classical results on Laplace 
and transmission problems in [§]). 

4.2 The Nernst Planck Equation 

Let us define V v , = F7 1 • V and ip e (., t) = ip e (X e {., 0, t), t). We consider the 
problem P2 in its Lagrangian version: For each i £ {1, . . . , J}, find 

A/" i>e (.,t) :=M,e(X e (.,0,t),i) 

satisfying 

dttiULt) = d 4 V Ve • fv Ve M, e + ^tii,eV Ve $e] e ) (£,*), (£,t) G n x (0,t ] 
M >e (e,o)=M,o,e(e), ?e«o, 

f^A4 )£ + ^M, e (£,t) = 0, K,t) G dOUdK X (0,t ], 

(4.76) 

where we have extended ip £ (.,t) to be zero in t G (to,? 1 ]- From the hypothesis 
H we have the compatibility condition 

d v Afi,oAO + -^vo,m(0 9»ihAZ> 0) — 0, £ G dO U ax . (4.77) 
Also, we consider the weak formulation for (|4.76l) : Find 

M, e £ L°°(0, t ; L 2 (n Q )) n L 2 (0, t ; ff 1 ^)) 

satisfying 

Y Sf^d^J (t) + d *J n J (v Ve M,6 + ^|M,eV Ve [^ £ ] e ^ • V v .cdedr+ 

/ / Mi^ t d£,dT = ( f Af i<0 ,e<;odz) 
Jo Jn \J£i J 

(4.78) 
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for all ? G H 1 ^ x (0,t )) and a. e. t G (0,t ). 

A direct calculation using (|4.68|) gives us that the problem (|4.76j) is uniformly 
parabolic. Using the regularity of the coefficients we can obtain the existence 
of a unique nonnegative solution G C 2+Q ' 1 (Jh x [0, t }) of (|4T76l - (|4"T7l) (see 
Theorem 5.3, Chapter IV in [25]) that corresponds to the unique solution of 
(|4.78[) . For sufficiently small to, we can obtain uniform estimates on A/i je . In 
fact, we have the following Lemma. 

Lemma 4.3. The solution Mi^ of \4.76\) - $- 77| ) satisfies 

sup f \NULt)\ 2 dt,+ (° ( |V ? M, £ | 2 ^dr<L 2 , (4.79) 
te[o,t ] Jn Jo Jcta 

where Li is given in |7~55p. 

Proof. First, we observe that a standard calculation gives us that, for all t G 
[0,*o] 

\ f (|M, £ (e,t)| 2 -|M, £ (C,o)| 2 )de + d 4 f f |V V£ M, e | 2 d^r = 
z Jfio Jo Jn a 



diZie '" 



< 



k bV Jo Jn„ 
d rt 



(4.80) 







V v M 2 dtidT + ^^ f f |M. e | 2 |V Ve [^] e | 2 ^r. 

S JO Jfio 



Now, from (|4.68l) . (|4.6ip and recalling the Remark 1, we have 



v Ve [^r(^)i 4 d£<i6 / |v c [^] e (e,oi 4 ^ 

4 



1G 



1G 



V e / ^(^r) 5e (t-r)dr 
o 



to 



V^ 6 (^,r) 5e (t-r)dr 



^ (4.81) 

de = ie/ |[v^ e (e,t)] e | 4 ^. 



From fil)?) and the definition of X e , we have ||F e (X e (., r, 0), 0, t)|| ,oo,O < 4/3, 
for each t 6 [0,to]- Then, using (|4.62|) and (|4.73l) . we have 

V |V Ve [^] e (e,r)| 4 ^dr< 16 / / t0 |[V^ e (C,r)] e | 4 dedT 
o Jn Jn Jo 

: !G / ° / |V c V £ (C,t)| 4 ^7 



JO: 

to 



1G 



< 



4096 



<*o 



f° f |F e (X E (x,r,0),0,r)V x V' e (x,r)| 4 dxdT 
Jo JS1,, T 

f° f |V x ^(x,r)| 4 dxdr 

JO Jf2 £T 

^(Li + \\ P o\\h, Ko + \\nh,oo) 2 - 



81 



(4.82) 
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Using the Holder's inequality and (|4.82p we have 

/ / |M, e | 2 |v V£ [^H 2 ^T 

Jo Jn 



<tl /2 ^{Ll + \\p4l 2 , Ko + \\nl 2 ,eo)x 



(4.83) 



\[° [ |V^ e | 2 ^ + ess sup \Wi, e {.,t)\\l, % n a 
I Jo Jsn te(o,t ) 



where for J\f i e we have used a well known Sobolev inequality (see inequality 
(3.8) in Chapter II of [25]). From (jTBgj) . (|4T53"|) and (gUt , w e can write P~gQ")) 



I \KA^)\ 2 dZ+ I I \S7^ he \ 2 d^dT 
Jn Jo Jn 

KT^aIT ( |V ? aL| 2 ^t + sup \\^A- Mla.no) + ^ — I l^ol 2 ^, 
yJo Jn a te[0M)] J * a min J n 

(4.84) 

for all t £ [0, to], where A and d m i n were defined in (|4.56[) . Recalling (|4.55[) and 
(H3S1), we obtain $£7^. □ 

Lemma 4.4. 77ie solution Af^ of J^.76'| )- f7- 77p satisfies 

f° I (d t U^) 2 dCdr + ess sup / |VM, e (£, i)| 2 d£ < C e (l + ||M,o||g,2,n ). 
Jo Jo te(o,t )Jn 

(4.85) 

w/iere C e = C e (£>, if , «i, «2, -Zi, ■ • . , Zj,di, e, K B ,0,^,t o , Li, L 2 ). 

Proof. Taking the L 2 -internal product of (|4.76|) with dtJ\fi, e , we obtain, after 
integrating by parts and the use of the boundary conditions, 

° (4.86) 

I -^ v ^ e • Vv.(ak#i, e )de) (*) = 0, 



/no 



for each i s (0, to]- Now, integrating by parts again, we have 

/ M. e V Ve [^ e ] £ • V Vc (d t M he )dA (*) = ( [ 9tK,e K, e V v AA} € • uds+ 

Jn J \JdOUdKa 

- / ft^i, e v v . • (a^Vv^k) (t) = f i / a t (M, e ) 2 v V5 [^] e • 

Jflo / \ Z JdOUdKg 

9 t M, 6 v Ve • (M, e Vv.[&] e )d£ ) (t) = (i / a t ((M, £ ) 2 v V£ [^ £ ] e • ^)d s + 



4 / (M, e ) 2 5 t (v Ve [v e ] e • - / a t M, £ v V£ ■ CA/^Vvj^n^ (*)• 

2 JdOUdKo Jn J 



(4.87) 



21 



Also, it is easy to check that 



Vv.AQ, e ■ Vv.fifrA£> = gd^VM.eHi - V Ve M,e • ((Fr^tVM.e). (4.88) 
Then, from (|4T87|) . fljjSp and integrating ([4T86]) in (0,t), we obtain 



(d t Afi,e) 2 dZdi 



o Jn 
2 



S2„ 



iF^Vf-A/^r^ ) (t) - y { j a K^MefdZ ) (°) 



Jo Jn 

((A4,i) 2 VvJ'0 £ ] e • vds) (t)+ 



2k b 6> 

2k b 6» 

dicZi 



28kb Jq 
d;eZ, rt 



dOUdKo 



dOUdKa 



((K,i) 2 V Ve [A} e ■ ")ds (0)+ 



(M l£ ) 2 St(V Ve [V e ] e -^)^dr+ 
9tM, e V Ve • (M, e V Ve [^] £ Rdr = 0, 

O 



for each t G (0, to]. Let us now obtain bounds on each term given above. 
First we observe that, from the incompressibility condition of v e , F^ 1 = (ay) 
corresponds to the adjoint matrix of F e , i.e., for each i, j s {1,2,3}, 



Oij = (-1) 



i+3 



V 9£ r 9£/ 9£ r 9£; 



(4.89) 



where X e — (Xi_ e , X2, e , X^). Furthermore, n, I are the largest elements of 
{1,2,3} satisfying n ^ i, I ^ j, respectively; m, r 6 {1,2,3}, m ^ i,n and 
r ^ Now, from (|4.63[) . 

X i , e (£,0,t) = ( i + f v ije (X e ((,0,r),T)d T , 
Jo 

for all (£,t) 6 O x [0, to] and i G {1,2,3}, where v e = («i, e , «2,e, «3,e)' As a 
consequence, if a;, = Xj )e , 



9X 



9& 



■(.,0,t) 



<8i 



0,oo,O 



t 3 

E 

fc=i 



9t>i, £ (.,r) 



0,oo,O 
3 



5X i>e (.,0,r) 



<1+ / ||V x v e (.,r)||o,oo,o V 
70 k=i 



96 

dX M (.,0,r) 



fir 



96 



0,oo,O 

dr, 

0,oo,O 



for all i,j £ {1,2,3} and t £ [0,t ]. Then 

ax 



E 



96 



-0,o,t) 



0,00,0 



< 1 + 3 / ||V x v e (.,r)|| ,oo,oE 



9X ij£ (.,0,T) 



96 



0,oo,O 
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which implies that 

dX„ 



-(.,o,.) 



0,oo,Ox(0,t ) 



<l + 3 / ||VxV e (.,r)||o.oo.oexp ( [ || V x v e (., r) \\o,oo.odr) dr 

Jo \Jr J (4.90) 

<l + C e J^ ||V x v e (.,T)|| ,2,oexp (°^J o HVxVe(.,r)|| 0!2iO dr^dT 



< C e {to, Li), 
using the Gronwall's lemma. Analogously, 
d 2 X, 



-(■,o,.) 



< Ce(to, Li), 



0,oo,Ox(0,t ) 



for all i,j £ {1,2,3}. Hence, 

||(5 4 F7 1 )(.,0,.)|| ,oo,ox(o,to) <C e (t 0) Li). 
Now we observe that, from (14.61)) and (|4.75[) . 

|V v .$2 >e ] e (£,t)| < / O | 5£ (t-T)||F e - 1 (f,0,t)V^ 2 , e (e,r)|dT 
Jo 

rto 

< C e / |V x ^2,e(x, r)\dr < C e (i )ess sup ||^2,e(-, OIIchoTT) 

Jo *e(o,t ) 

< C 6 (ki, k 2 ,O,K ,j, Z r , . . . , Zj,t ,L 2 ){l + ||po||o,2,jr + ll*||i,2,ao), 
for all £ £ H,q and a. e. t £ (0,to). Similarly, we have 

|V Ve -(V v J^] e )(e,r)| < / Q |.g e (t-T)||A x ^,e(x,r)|dr 

Jo 

< C e (K 1 ,K 2 ,0, Kq,J, Z 1 ,...,Zj,t ,L 2 ){l + ||/0o||o,2,Ko + ||*||l,2,eo), 

and 

|ft(V Ve [^] £ )(e,t)| < [° \d t9e (t-T)\\V x TP 2 , e (x,T)\dT 
JO 

<C e (t Q )ess sup ||^2,e(-)*)|lci(n~7) 

t€(0,t o ) 

< C e (K 1 ,K 2 ,0,K Q ,j,Z 1 ,. .. ,Zj,t ,L 2 )(l + ||po||o,2,if + ||*||i,2,eo), 

for all £ £ fio and a. e. t £ (0,to)- 

Using tgZBBf, P~§H and dUTSJ), we obtain the estimate 

f f V Ve M, e • ((F~ 1 )iVfA/i |6 )d$dT 
Jo Ja 

<C e (to,Li) / ||V s M, e ||^ rfT<C e (to, J Li,£ 2 ), 
Jo 

for each t g (0, to]. If 1/2 < s < 1, the trace theorem and (I4.92j) give us, 

dOUdKa J 



(4.91) 



(4.92) 



(4.93) 



(4.94) 



< a 



00 



(N^fds+ / {N^fds) <C e \WU-Ml%a , 



dK 
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where C e = C e (/«i, k 2 , O, Kq, e,7, Z\, . . . , Zj,to, po, Zr 2 ). Hence, from inter- 
polation results (see Theorem 1.3.7 in [31]), Young's inequality and (|4.79p . we 
have, 



((M, £ ) 2 v v jv e ] e ^)ds (t) 

<c e |i^ iie (.,*)ii^iiM, 6 (.,t)ii?r 2)no 

<C e + S p 1/s ||V e M,e(.,*)||^ 2 ,n , 

where p > and C £ = C e (e, s, p, k 1? k 2 , C, #o,7> • • • > Zj, t , p , L 2 ). Let 
us to recall that suppWi.o.e CC Qq, so that 



((M, e ) 2 V Ve [^] 6 -i/)ds (0) = 0. 

dOUdKa J 

Also, using (|4.94|1 . trace theorem and (|4.79|l . for each t £ (0, to], 

f / (M, £ ) 2 a t (v Ve [^] e • v)dsdr < c e [ \\MA-,T)\\l 2Mo dr < c e , 

JO JdOUdKo Jo 

where C c = C e (e, n x , k 2 , O, Kq, 7, Zi, . . . , Zj, t , po, £2)- 

Furthermore, using the Young's inequality, (|4.92[) . ()4.93p and (|4.68jl . 



/ / 9*M, e V Ve •(M >e V Vt [-0 e ] £ )^dr 
Jo Jq 

< \ I I (dtK.e) 2 d^dr + l f f |V Ve -(M, e V v J^ £ ] e )| 2 d£dT 
^ jo Jn * Jo Ja 



<\ I I {dtMi,e) 2 d£dT+ I I |V V£ M,e| 2 |V Ve [^] e | 2 ^dr- 
^ Jo Jq Jo Jo,q 



- f f |AU 2 |V Ve ■ (V v AAY)\ 2 dtidT 
Jo Jn n 



Jo Jn 



2 Jo 

< 



1 '* 



2 



n Jo Jn Jo Jn 

2 . 



{d t Afi,efd^dT + C e , 

Jn 



where C e = C e (e, Ki,k 2 , O, Kq,^, Zi, . . . , Zj,to,po,~^,L2). From the previous 
results, choosing p = ( g^f^- J we obtain 

i f / (d t Af i , e ) 2 d£dT + ^ [ |v c M, £ (t,OI 2 ^<c e (||M,ollo,2.n„ + i) 
* jo Jn 10 J(] 

for each t £ (0, to], which implies in (|4.85p . □ 

4.3 The Navier-Stokes Equations 

In order to deal with the problem P3, we define 

J 

F e (.,t) = -e^(^/Vi, e Vv>e)(.,t) 



i=l 
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and consider the linear problem related to (|3.34p in Lagrangian coordinates: 
Findu e (.,i) := u e (X e (.,0,t),t) such that 

n(^ Pi0 u e • d t ip - 2nD v ,(u e ) : D vC (ip)) d£ dr+ 
J 

+ / V P ,oUo,e ■ <f(0) d£= I Mp, u e ■ ip dn (t) + / / -p f ¥ e ■ tp d£ dr, 

JO \Jo J JO Jflo 

(4.95) 

for a. e. t £ (0,to) and V</? € S e , where 

S e = {ip e H\0 x (0, t )) 3 , V v « • <p = 0, <p\ go = 0, 1/p.oAr-V = in x (0, i )}. 

Lemma 4.5. TTie problem has a unique solution 

u e € L°°{0, t ; L 2 {0)f n £ 2 (0, to; Jtf(O)) 3 , 

swc/i t/iat V v e • u e = 0, /^ Pi o^v»(u e ) = a. e. in O x (0, to). Furthermore, u e 
satisfies \3.3J$ and we have the estimates 

ess sup ||u e (., t) || 0,2,0 "I" 

/ llV^llg, 2>0 rfr < Lf, (4.96) 
te(o,t ) Jo 

where Li is given in |^.5ff| ), and 

/•to 

ess sup ||u e (.,t)||? )2j0o + / ||^Ue||o,2,ei dr 

te(o,t ) Jo ^ g7 ^ 

+ / ° K||L,n f A < C e (||uo||g AO + Ll + ||po||2 A jr + ll*lll,2,o) , 
Jo 

where C e = C e (r),Ji p ,~pf ) K , <D, Zi, . . . , Zj, e, Ci, C**). 

Proof. The existence and uniqueness is based on a special Galerkin approxima- 
tion technique and was established in [TS] jointly with the estimate (|3.29|) for 
u £ . Now, from Schwarz inequality, Sobolev embedding (see (|4.57p ). (|4.73l) and 
(|4J9]), 

f° [ |F e | 2 dxaY < 2e 2 f° [ ^ \Zi\ 2 \Afi^\ 2 \Vip e \ 2 dxdT 
Jo Jn CtT Jo Jn tiT i=1 

< 2e 2 [ ° ]T l^| 2 HM, e |lo.4,a l! V V^llo.4,a e T *r 

<2e 2 Y,\Z*\ 2 Ctess sup ||V^ e (.,t)||o,4,n. T f " ll-AMi,2,QcA 
i=1 te(o,t ) ' Jo 

J 

< 2e 2 Y,\Z l \ 2 ClL 2 C 2 (L 2 + \\po\\l*K + ll*ll? A ao)- 

»=1 

Then, recalling (ET2"9"1) , the estimates (|4~M1) and (|4~551) give us QjZMty . 
From the regularity results established in [15], we have 

rto 

ess sup ||u e (.,t)||f j2 ,no + / \\ d t u e\\o.2,O dT 
te(o,t ) Jo 

+ j£ ° KHI^aA < C« (lMlo,a,o + j[ " l|F e ||L,o 5 , T rfr 



(4.98) 
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where C e = C e (r],]I p ,~p f , K , O). Hence, (|43T| follows from □ 
4.4 Fixed Point Procedure 

Let us recall that we have chosen (v, . . . , G £> x <Y x . . . x X. Now we 
observe that 

u e (x, t) = v C:f (t) + w e (i) x (x - x £ (£)) 

for (x, i) € K ty t x (0,to), x £ (i) is the center of mass of K e .t- Let us define the 
map 

g £ (v, = (^^.ej 1 ^),^,,,...,^), 

where (x £ ,0 £ , 6^" x (y 6 )) G F to , x e (t) = x c (0) + /„* v Cie (r)dr, £ (i) - /* w e (r)dr 
and 6^ 1 (y e ) denotes the incompressible component in flo of 6 _1 (y e ); as de- 
scribed in [15] . y £ is a divergence- free vector field constructed from u £ , v £ and 
from suitable rigid current functions. In order to show that this map has a fixed 
point in B x X x . . . x X we need to prove some technical results. 

Lemma 4.6. G e is a continuous map from B x X x . . . x X into itself. 

Proof. First we observe that G e applies B x X x . . . x X into itself. In fact, using 
(|4.79[) we have the estimate for A/i i£ . As in [T5J, the estimate for (x e , 0^" 1 (y c )) 
follows from (|4.96l) . Now, let us to consider (v, <&t, . . . , dj) E B x X x . . . x X 
and sequences v„ — > v in B, di. n — > di in X, as n —> oo. We need to show that 

G £ (v n ,#i,„,.. . ,f?j, n ) -> G £ (v,#i, . . . ,tfj) in B x X x ... x 

as n —> oo. As before, we denote v n = 0(v n ), v = 8(v), v„ )£ = 1Z e {v n ) and 
v e = 1Z e (v). Also X„ )£ and A £ represent the Lagrangian flows of v„ j£ and v £ , 
respectively. From the properties of the function (see [H]), for fixed e > 0, 
v„, £ -)• v £ in L°°(0,to;L 2 (0)) 3 D L 2 (0,t Q ; H^(0)) 3 , as n -)■ oo. Henceforth, 
denoting 

F n>£ = VX e ,„ and F £ = VX £ , 

we have 

||(F- 1 e -F e - 1 )(.,0,.)|| , oo , Ox(0 , t0 ) 

< C(e,t , £>,Li)||v„ )£ - v e || L oo (0!to;i 2 (c ,))3 ni 2 (0! i . H i( C ,))3 0, n ->• oo. 

(4.99) 

In fact, this follows from (I4.63I) combined with (I4.89|) and (|4.90|) . 

Now, we define 

(0<,»,e(->*) := ^,„,e(X„ )e (.,t,0),t),^, e (,t) := 5 i , e (X e (.,i,0),t)), (4.100) 

and consider the Lagrangian versions of (|2.3p - (|2.9p . for a. e. i 6 (0, to), 

J 

V Vn , e -(k(.,0)V v „,> n ,e) = -47re^g iin)e -47rpe in (H l {0))\ 

i=i 

^n,2,e(£,*) = V>n,l,e(£,t), £_ E dK , 
(Kid v 1pn,l,e - «2#i^n,2,e)(£,*) = 0, £ G <9if , 
&,2,e&*) = *«(£)> 

(4.101) 



2G 



V Ve .(k(.,0)V v .&) =-47re^§-, e -47rp e in {H l {0)) f , 

t)=&,e(£,t), «.e.^aif , (4.102) 
(/ti<9„^i, e - K 2 d u ip2,e)(^, t) = 0, a. e. £ 6 9iT , 



where, 



and 



P0,e(O> £ G ^0, 

o, e g no. 



Considering f n<e := ^ nie — ip e as a, test function in (|4. 1 1[) and in (|4.102l) . 
after some calculation, we obtain 



«o|V VBi ./n, e rdeJ (*) < (J «o|(Vv„.> £ - V v > £ ) • V Vn ,./«,e|d£ ) (*)+ 



+ K \(W Vn J n , e ~ V Ve /„, e ) • V Ve ^e|^) (*)" 

+ 47r e J] |Zi| Qf |^,„,e-^,eli/n,eMe) (*) 



< IIO^-lT^O.OIIo.^Oxto.to) 

+ ||V € / r , !e (.,t)||o,2,o||V Ve -0(.,t)||o,2,o) + C||(^^, e - ^, e )(.,t)||0,2,o||/^, e (-,*)||0 !2 ,O, 



where C — C(e, Zi, . . . , Zj). Clearly, as ||($i, n — *)||o,2,o — > with n — > oo, 
we have 



|V 4 / n , e (.,i)|X (t) -)■ 0, n -»• oo, a. e. * G (0,i ) (4.103) 



using (|4.73l) . Poincare's inequality, (I4.68|) and (|4.99|) . 

We define (Mi, n , e ,^fi,e) an( l (u„, e ,u e ) as the solutions of (|4.78|) and (|3.34[) 
considering ([^„, e ] e , v„, e ) and ([-0 e ] e ,v e ), respectively. Taking W iin , e := M,n,e - 
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Ni, t as a test function in (j4.78[) we obtain, after some calculation, 

/ (Ui, n , e (£,t)fd£ + di f f |V v „ E W w | 2 d^- 
Jn a Jo Jn 

= -di [ f (V v „, e A/' i , e - V V€ Af i>e ) ■ V Vn ,Mi,n,ed£dr+ 
Jo Jn B 

-di [ [ (V Vn< Mi,n,e ~ VvM,n,e) ' V Ve M,e^T+ 

Jo Jn 

f f ^,„. e V v „_ e [^„, e ] e • V Vnt Mi,n,ed£dT+ (4.104) 
o Jn 



diZie 
k b 



diZie '" 



43- w 



M, e [V v „ JV™,e] £ " V v „ e [V e ] £ ] ■ V v „ e U l .n,ed£,dT+ 

n 



kbO Jo 

d,Z, r >< 



f f M, e [V v „, e [ipeY ~ V Ve [^ £ ] £ ] • V Vn> Mi,n,ed&T+ 
JO Jfln 



kbO Jo 



Q 

n 



for each t € (0, to)- Then, using Holder's inequality, Sobolev embedding and an 
analogous estimate as in (|4.83p and (I4.84[) we have 



/ {Ui, n , e (£,t)) 2 d£+ f f \W C Ui tn , e \ 2 d^dT 

Jn Jo Jn 

< CHOP- 1 , -F7 1 )(.,0 J .)|| 0( oo ) Ox(0,to) x 



||V S A/" i ,e(- ) T)||o,2,nol|V e Z4,n,e(-,T)||o,2,n dT+ 

+ / ||V f [^]%, r)[| ,4,fi IIM.eC-, -r)|[i,2,o ||V £ W i ,„, e (. ! r)|| ,2,a o dr 



+ 



0,oo,n x(0,t ) / ll^cdV'n.e] 6 — [i>e] 6 ) (■ , t) \\ 0,2,0 \ | VgWi,„, e (■ , r) || o,2,n dr, 

Jo 

(4.105) 

where C = C*(O , k b ,0, Z u e, di). Hence, using (j499)) . (1462)) . (|4~T03)) and (|4J9)l 
we obtain the convergence 

ess sup f (U l , n ,e{Lt)) 2 d^+ f f |V^,„ :e | 2 d^r^0, n^oo. (4.106) 
te(o,to)Jti Jo Jn 

Now, we define 

,/ J 
F e (,t) = -e 5^(2^6 Vv.Ve) (-,*), Fn,e(-,i) = -e$3(^, n , e V Vn >n, e )(.,t), 



1=1 1=1 



respectively. Taking q„. £ := u„ !£ — u e as a test function in p. 341) . using (|3.29[) 
and we have q n , e -> in L°°(0, t ; L 2 (0)f n L 2 (0, i ; ^oHC)) 3 if w e can 

show that ^ 

/ ||(F^-F £ )(.,r)|| 2 . 2iao dr^0 
Jo 

with n — ?> oo, for a. e. ( £ (0,to)- This result follows writing the above 
expression as in the right side of (|4. 104|) and from the estimates obtained in 
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gTD5l) . Then, following [15], 

(v,^ £I e 2 1 (y«))^(x e J £ ,e 2 1 (y £ )) m 6, 

as ?i — > oo. Hence, recalling ()4.106|) we obtain the continuity of G c . □ 

Proposition 4.1. The operator G e has a fixed point in B x X x . . . x X . 

Proof. First we observe that G e is a compact operator. In fact, this follows from 
Lemma 14.61 the rigidity properties in Kq and from the Sobolev embeddings in 
fio for Mi, c (as a consequence of the estimates (|4.79p and (|4.85p ) and u e (from 
the estimates (|4.96p and (|4.97[) ). Then, the map G e has a fixed point in B, using 
Lemma S3] and Schauder's theorem. □ 



The first time step on (0,t ) is completed. From interpolation results (see 
Theorem 1.3.8 in [3T]), u e (.,to) € i? 1 (C) 3 . Therefore, we can proceed similarly 
on the interval (to,ti), considering (u e (., to), A/" e (., to)) as the initial conditions 
for the problems (|4.95[) and (|4.78p . respectively. We give below some details 
related to the obtention of the bound (|4.79[) in this second step of the procedure. 

Let us consider ip e and A/i, e the solutions of the problems PI and P2 in 
(0,to), respectively. We extend these functions in (to,ti) to be the solutions 
of PI and P2 in this interval. First, we observe that (|4.73p is valid in (£o,£i) 
using the Lemma 14.11 From a similar argument as in (|4.80p . we have, for all 
t 6 [t ,*i], 



\ [ \MUU)\ 2 d^ + ^ f [ \V^\ 2 dtdT 
<JS/V IM,e| 2 |V Ve [^] e | 2 ^r+ (4.107) 



'to JQq 

\K(tto)\ 2 d£ 



and, for all t G [0, to], 



1 /' • ... . 

<SS/ / |M,e| 2 |V Vt [^] e | 2 d^+ (4.108) 



/ \JViAS,t)\ 2 dt + % [ f \VzM, e \ 2 d(dT 

JQo ° JO iSiri 



2k%0 2 jo Jn 



\^,0)\ 2 dt 
n 



Considering t = t in (|4.108p . we obtain the validity of inequality (|4.108p for all 
t £ [0,£i]. As a consequence, we obtain a version of (|4.79l) in the interval (0,£i), 
following the lines of the proof of Lemma l4~3l In particular, f|4. 79[) is valid if we 
consider the interval (to,t\). A similar argument can be used in order to obtain 
the estimate (|436T) in (0,£i). 

Repeating all the process in each step (ti-i,ti) (recall that the number of 
steps N = T/to only depends on e and L\) we obtain, for each e > 0, (u e ,M e , tp e ) 
as the solutions, respectively, of the problems (|4.95[) . (|4.78[) and (|3.3T[) in (0,T) 
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and /i c that satisfies (|4.71[) . If we consider u e and Afi jC in terms of the Eulerian 
coordinates, we see that, for each < e < e', (u e , fi e , A/i, e , . . . , A/j )e , ip e ) is a 
solution of the problems PI, P2 and P3 and we can obtain an approximate 
sequence (u„, A/i, n , . . . ,Afj, n , ipn)neH of solutions for (|3.34|l - (|3.38|l . Using 
(fOSj) . P~79")) and (|4T55|| in each time step (tj-i,**) we obtain, for all n G N, 
i S {!,..., J}, 



(4.109) 



eSS Sup ||M,n(-,*)||o,2,Q„ t + / II VM,n||o,2,Q„ A < 4i 2> 
t6(0,T) ' JO 

ess sup ||u n (.,f)[[g 2) o+ / ll Vu ™llo,2,G^ T < L i- 
te(o.T) Jo 



«6(0,T) 

Also, 

ess sup ||^„||o.oc.,D < maxj/Ij,^}, 
te(o,T) 



(4.110) 



ess sup max(||VV'i,„(.,i)||o.4,/f„ t , ||V-02,n(-,Ollo,4,o„ t ) < B*, 
*e(o,T) 

where B, = C 2 1/2 (i2 + ||p ||g, 2) x + ll*HUao) 1/2 and 

inf inf dist(if nt , 90) > 7 > 0. (4.111) 

nSNte(0,T) ' 

Remark 4.2. In order to establish convergence properties for Af n it is conve- 
nient to consider its extension to fixed domains. So we consider the standard 
extension operator E : ff^fio) — > i? 1 (M 3 ) such that Ef = f a. e. in fio, 
HV(J5/)||§ i3>R a < C||V/||g i2ino , \\Ef\\l < C\\f\\l 2A) and £/ > a. e. in 
R 3 as / > a. e. in Hq; the constant C depends only on 8Kq and dO (see 
Chapter 5 in [17]) • Then, from (|4.109j) and (l4~68|) . we have 

ess sup \\K: n (-M.2,R3+ C \\VKn\\l.2^dT<ACLl (4.112) 
te(o,T) Jo 

where A/> n (x, t) = AA*„(A„(x, t, 0), t) and A^*„ = 



5 Convergence Results 

As a consequence of the bounds (|4.109j) . (I4.110[) and (I4.112p . there exist subse- 
quences {u nfc } fceN , {n nk }km and {A/> nfc } feeN such that 

u„ fc weaklyin L 2 (0, T; ^(O)) 3 , 

Hn h ->-/i weakly* in L°°((0, T) x £>), (5.113) 
^*„ fc ->A/"* weakly in L 2 (0,T; 

for some (u, n,Af*, ■ ■ ■ From the compactness results for linear transport 

equations of Di Perna-Lions (see [IB]), there exists a subsequence {/Un*, . }jeN 
such that 

\i nh . —> jj, strongly in C([0, T};L P (D)), (5.114) 
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for all 1 < p < +00. Moreover, from (|4.109|) we can show that u e L°° (0, T; L 2 (0)) 3 
and N* £ L°°(0, T; L 2 (0)), using a uniqueness argument. For simplicity we de- 
note (u„,jU„,A/£ n ) instead (u„ fc _. , fj, nk . ,jV£ nk . )■ 

It is a routine to check that (/i, u) satisfies (|3.35j) and (|3.36[) and u satisfies 
(|4.109j) . In particular, from (|4.111j) , this implies that (|2~2]) is valid for the 
domains K t and fit that correspond to u. Below we establish some compactness 
results. 

Lemma 5.1. For all h > sufficiently small, 

sup/ / fi n \u n (x,t + h) - u„(x,i)| 2 cftdx < C7i 2/5 , (5.115) 

nGNJO JO 

w/iere C = C(T, L 2 , L 1; 6*, K , e, O, Zi, . . . , Z n ) and we have extended u„(., r), /!„(., r ) 
<o &e zero if t > T . 

Proof. Let us to denote w n (.,t, h) = u n (.,t + h) — u n (.,t) for all t > 0. As a 
consequence of (|4.11ip and considering the construction and notation introduced 
in [15], it is possible to find ±I^(w n (., r)) £ S n (r), Vt £ [t,t + h], such that 



||w n -n^(w„)||g iPj0 dt < C p f3 



l-s 



l Vw n|lo,2,0^> 



/ ||vn^(w„)||2 pi0 dt<cr s / ||Vw n ||g )2>o( ft, 
Jo Jo 



(5.116) 



for s = 3(1/2- 1/p) ifp£ [2,oo), s = 3/2 if p = 00. 

If we extend Zq iit (.,t), A/i,„(.,r), -0n(-,' r ) and II^(w n )(.,T) to be zero if 
r > T, we have 

f T f f t+h 

/ / / ln n/r F n (x,T) -ILp(w n )(x,t)dTdxdt 
Jo Jo Jt 

J pT pt+h p 

<C(e,0)y2\Zi\ / ||Vn^(w„)(.,t)|| .oo.O / {\ln n Mi,nWn\){*,T)dTdxdt 

i=1 Jo Jt Jo 

<C(e,0)J2\Zi\( [ 
i=l \ J o 



\V2 

l|vn^(wj|| 2 iCOi0 dt x 



t+h 



{\In nT Afi,n\\Vip n \)(x,T)dTdx 



1/2 



dt 



j / pT r ^-t+zi ] r pt+h ] \ 

^C^CLOr^El^l (/ I IK«llo A n B<T dr jf ||V^]]^ 2)0 dr dil 
< C(e, O, Li, Zi, . . . , Z 7 , L 2 , B„ T)^' 2 h < Ch 2 / 5 

choosing f3 = ft 2 / 5 ; where we have used the Schwarz's and Poincare's inequali- 
ties, (|5.116p . (|4.109[) and (|4. 1 10[) . As a consequence, following the steps of the 
proof of Lemma 4.1 in the reference [15] (see also [13]) we obtain the inequality 
(15~TT5)) . □ 



1/2 



31 



Lemma 5.2. For all h > sufficiently small and each i G {1, . . . , J}, 

sup/ / ln nt \Af* n (x,t + h) -N* n {x 7 t)\ 2 dxdt < Ch 1 / 8 , (5.117) 
new Jo Jr 3 

where C = C(di, . . . , d n , L2, L±, C, T, kb, 0, e, Z\, . . . , Zj, O, B*) and we have 
defined Iq u t (.,r), N* n {.,r) to be zero if r > T '. 

Proof. Following [13] , we have a scalar version of (|5.116|) . More precisely, (|5.116|) 
is valid for functions in H l (M?). In this case, ILp(f) = f*T]p, where r\ G Cq°(R 3 ) 
satisfies J R3 rjdx = 1 and, for each /3 > 0, rjp(x) = /3~ 2 ?7( / 3 _1 x). 
Let us consider i G {1, . . . , J} fixed and denote 

w*,„(,*,/i)=A/;y,* + /i)-^y,t). 

Then, for each 11 £ N, we bound by 

/ / Xn nt (x)W? n (x,t,/i) ( ix^<P 1 +P 2 +P3 ! 

JO JR 3 

where, 

Pi ■■= [ T [ (2b»,*+ fc -Xo„, t )(x)(^ n ) 2 (x,i)dxdi, 

JO JR 3 

P 2 :=/ / (^„ it -Tn n , i+ J(x)(/V;%) 2 (x,t + ft)dx^ 

JO JK 3 

^3 := [ T f (la n , t+h (x)A^*„(x, t + h) - ln n , t+h (x)A/^(x, t))W i ,„(x, i, fc))dxdt. 

JO JK 3 

Using (|4.112p and Lemma 1 of the reference [13] , we have 



Pi < C/3 1/4 + C 



(2n„. t+fc - 2n„, t ) (x) (IL3 (A/^)) 2 (x, t)dxdt 



(5.118) 



where C = C(L,2,C,T). Now, extending lZ n (u n ) to be zero in M 3 \0, we see 
that T n (x, i) = Xn„ t (x) is the solution of the transport problem 

d t In + V • (l n TZ n {n n )) = 0, !„(., 0) = 2n . 

in P'((0,T) x M 3 ). Then, using (|5.118j> . we have 



T r r t+h 

V • {2 n ll n (u n ))(x,T)(Il (N* n )) 2 (x,t)dTdxdt 

Jk 3 Jt 



Pi < cp l/i + c 

pT pt + h p 

+ C / |(X n 7^(u n ))(x ) r)||(n^(A^* ))(x,t)||V(n^(^* ))(x,t)\dxdTdt. 



Now, from Holder's inequality, the Sobolev embedding theorem, (|4.64[) . (|4. 109[) . 
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(|5.116p and 1)4.1120 . we have 







| (Z n K n (u„)) (x, r) 1 1 (IL5 (A^J) (x, 1 1 V(n^ (A? J)(x, i) |dxdrdi 



/ | (Z„ft„ (u„ ) ) (x, r ) 1 1 (H^ (M* n ) - M* n ) (x, * ) 1 1 V (11^ (M* n ) ) (x, t) I d*drdt+ 
/ | ( l n K n (u n )) (x, r) | |JV? n (x, t) 1 1 V(n^ {M* n )) (x, t) | dxdrcft 

< / T r +/l ||^„(u n )(.,r)|| O!6iR3 ||V(n (A4:„))(.,t)|| o ,3, R 3X 
Jo Jt 

x ||(n^(JS^ n )-^ n )(.,i)|| 0l3 ^sdrdt+ 

\\TZ n (u n )(., r)||o, 2 , K 3 1| V(n (3 (A4*„))(., t)|| ,3,K 3 ll^n(-. *)llo,6,R»drtft 

< C/i 1 / 2 ^ ||u„(.,r)|| 2 i2 ^^ (j* ||V(n^(^J)(.,t)||2 3)R 3^ x 

x ^IKn^^j-MtJl.^)!! 2 ^^ + 

where C — C(L2, L\, T, C). As a consequence, 

Pi < C((3 1/4 + h 1 ' 2 ^^ + hp- 1 '*), (5-119) 

A similar estimate can be obtained for P 2 . Now, an elementary calculation show 
us that 

P3 < 

If ((I„A^* n )(x, t + h)- (X„A^* n )(x, t)){W i>n - n^(W il „))(x, t, h)dxdt 

T 

((2„^ w )(x,t + ft) - (2 n ^ w )(x,t))n^(Wi,»)(x,t s /i)dxdt 

Denoting the first an the second terms in the right side above by P3 and Pg', 
respectively, we have 

1/2 



( ,t X 1 / 2 (5.120) 

x n \\(i n M*j(.,t+h)-(i n K: n )(-Mh^ dt ) 

^(P^C)/? 1 / 2 , 

using Schwarz inequality and (|4.112p . In order to estimate Pi?, for t 6 [0,T] 
fixed, we take q(x, r) = n^(Wi in )(x, t), (x, r) G R 3 x [0,T] as a test function 
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in (|3.38p . Then, we have, 



+ 



(XJv^ n n /3 (W i ,„))(x,i)dx+ / / n j9 (W i ,„)(x,t) (7e„(u n )-VA^:„)(x,r)dxdT- 
3 Jo ./n n , T 



AT ii0)n (x)n^(W i>n )(x,t)d: 



x 



and 



(I„A^t„)(x, t + ^)n /3 (W i ,„)(x, i)dx+ 



o Jn„ 

t+h 



n^(Wi,„)(x,t) (7e„(u„) • VA^*„)(x,r)dxdr+ 

+ * i + X ( VA/ ^" + S^" V ^" ]n ) (X ' T) ' VII ^ Wi '«)( X ' *) dxdr 
M, ,„(x)n^(W 4 ,„)(x,t)dx, 



where we have extended [?/'«]"(■, t) to be zero if t > T. Subtracting the above 
expressions, we obtain, after some calculation, 



t+h 

d T {l n Af* n ){x, r) n^(Wi, n )(x, t)drdx+ 
t+h /■ 

/ n fl (Wi,„)(x,t) (ft„(u n ) • V/V£ n )(x,r)dxdT+ 

t ^f2 ra ,T 

/•t+A /■ 

di / / (V/V£ n + Af* n V[i> n ] n ) (x, r) • Vn^(Wi, n )(x, t)d X dr = 0. 



Hence, 



t+h 

d T {l n K*n ) (x, r) (Wi, n ) (x, t)drdxdt 



IQ JR-i Jt 
cT r t+h 



< / |n /3 (W i ,„)(x,t)||7e„(u„)(x,r)||V^„(x,r)|dxdr^+ 

J o Jn„, T 

/ |VA^ n (x,r)||Vn / s(W i , n )(x,t)|dxdTdt 
Jsi„, T 

+ dl2 Se Jo Jt +h L l-^t»( x ' r )ll v ^»] n ( x ' T )ll vn ^ w ^)( x ^)l dxdTdt - 



From Holder's and Poincare's inequalities, (|5.116[) . (|4.112|) . (|4.109|> and (|4T64|) 
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we obtain the estimates 

r T r t + h r 

\ |n^ {W i>n )(x,t)| |ft„ (u„) (x, r) 1 1 VA^t„ (x, r) | dxdrcit 

fT rt + h 

< / / ||n fl (W < , n )(.,t)||o, 00l R3||», l (u n )(. J T)||o,2,o||V^ n (. J T)||o I 2,odrdt 
Jo Jt 

T ( r t+h \ 1 / 2 / pt+h 

ix „ i ii,, m .. ; i f 1 ~ ,M " 

\Jt 

1/2/3-3/4 



<C(L 2 ,C)J[J^ |]7^(u n )(.,T)||^ dr] (/ ||V^„(.,r)||g )2iC ,dT 



(5.121) 



(5.122) 



rT rt+h r 

/ / / |VA^*„(x,r)||Vn /3 (W l ,„)(x,i)|dxdrdt 
Jo Jt Jn n ,T 

< £ U t+h \\^U x ' T )\\l2,odr) x 
t+h \ Va 

||vn j 9(Wi,»)(x > *)||g iai0 dT) df 

<C(L 2 ,C,T)J^ ^ + |lVn^(W i ,„)(x,t)||2 20 dT^ dt 

< C{L 2 ,C,T)h 1/2 , 

' f t+k f |A4: n (x,r)||V[^ n ] n (x,T)||Vn^(W l , n )(x^)|dxdrdt 

< J ||Vn^(W i ,„)(,t)|| ,oo, R 3 J \\K*n\\0,2,o\M^n} n \\0,2,cdTdt V ' ' 

< C{L 2 ,B*,(D)(3- 3/i h. 

Choosing (3 = h 1/2 and using (|5.119|) - (|5.123p we obtain the result. □ 

From Lemmas 15.11 and 15.21 and from Frechet-Kolmogorov Theorem (see The- 
orem IV. 25 in [6]), up to a extraction of a subsequence, 

u n -m strongly in L 2 (0 x (0, T)) 3 , 

V V " (5.124) 
A/? n -KA/? strongly in i 2 (0 x (0, T)), Vi G {1, . . . , J}. 

In particular this implies that N* > a. e. in R 3 , as it is easy to check. Now, let 
us consider the problem p.37[) with the data (K t ,Q, t ,Af*, . . . ,J\f},p, From 
Lemma |4~2"1 this problem has a unique solution ip(., t) G i? 1 (C) that satisfies the 
bound (|4.73p . As a consequence we can prove the following convergence result 
for ip n (.,t). 

Lemma 5.3. Let consider (if) n ,i/j) the solutions of the problem \3. 37\ ) corre- 
sponding to (K nit ,£l ni t,J\f^ n , . . . ,^Vj >n , Pn,*n) and {K u Vt u M*,...,J\fj,p,^), 
respectively. Then, 

[ [ |V(V>n - ^)| 2 dxdi 0, n-+oo. (5.125) 
Jo Jo 
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Proof. First we consider the standard extension operator 

E : H^dO) -> H\0), 

i.e., E (f)\dO = f and \\E (f)\\i,2,o < Co||/||i,2,ao. 

Then, defining $ = E ^ and $„ = E ^ n we have from (|4.65p . 

||*n||i,a,o < Co||*||i,a,oo. (5-126) 
Furthermore, it is clear that if T„ := ty n — ^, 

||T„||i, a ,o-»-0, ti y oo. (5.127) 
Note that tp n — ^„ + g L*; then, from (|3.37|) . we have 

(J KVtp ■ V(^ - 1pn ~ $ + $n)dxj (t) + 

- 47re^ Qf Z*A/?(^ - Vn - $ + $„)dx) (t)+ (5.128) 

- 4tt (7 - Vn - * + (t) = 0. 
Similarly, observing that ^ — $ + ^ n G , we have 

K„V^„ • V(V>n -i>~ *n + *)dxj + 

Zi r n (K*n) (?l>n-i>-§n + ^)d^J (t)+ (5.129) 

4^ ( / p„(^n - 1p - $n + *)dx J (t) - 
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Let us set A n , 4 = K n ,t D i^t, i?„, t = f2„,t n fi t and f n = ifr — ifrn- Then, summing 
(|5.128|) with (|5.129j) . we obtain, 



J^n n \Vf n \ 2 d^j (t) = 



(k - k„)V^ • V/„dxj (t) - Uf ( K - K „)W • VT„dxJ (t)+ 
jf K „V/„ • VT„dx^) (t) + 47re^ f jf Z,(AC - r n (A/£ n ))(/ n + T n )dx^ 

4^e V f / Z*A/?(/„ + T„)dx j (t)+ 

(t) + 4tt f / (p - Pn )(f n + T n )dx) (t). 



4 ^ e I] f / ZiT n {M* n )(f n + T„)dx 

j=1 \JK t \A n , t 



3G 



Holder's inequality gives us 

J^\Vf n \ 2 d^j (t) < 
< C(\ Kl - «2|||V^(. J t)||o,a,jc B . t \A B , t (||V/ n (. J t)||o J 2,o + l|VT„|| ,2,o)+ 

+ |«1 - «a|||V^(.,t)||o,2,Jf t \A l ,.(l|V/n(- > *)||o,2,0 + ||VT„|| ,2,o)+ 
,7 

+ l|VT n ||g j2i0 +^||(^-r B (^ n ))(.,t)||o,2 1 o||(/n + T n )(.,t)||o,2,0 



i=l 



+ E(ll r "(^t«)(-'*)ll0.2^A^, t + ll^(.,t)||0 J 2,K„ A A„J||(/n + T n )(.,t)|| ,2, O + 

4=1 

+ \\{ P - /3„)(.,t)|| ,2,A„, 4 ||(/« + T n )(.,*)||0,2,O + 

Q,3/2,K t \A nit )(ll(/n+T B )(.,t)|| 0) 3, O )), 

where C = G(Z\, . . . , Zj, e, fti, k-i). Using the Poincare's and Holder's inequal- 
ities, the Sobolev embedding, (|4.110l) and (15. 126|) we obtain 



J \Vf n \ 2 dxj (t) < C(B» + \\ni,2,do){\K n , t \An,t] y4 B* + |^\^n,i| 1/4 B*+ 

.7 J 

+Eii(^- r «W)(-»*)iioAo + Ei^AA., t r /4 ii^?(^)iUAo+ 

4=1 4=1 

+ X;i^t\A.,t| 1/4 |kn(A^ Tl )(.,t)|| ll2 , + 

4=1 

+ ||( / 9- / 9„)(.,i)llo,2^„, t + |^A^n,tr /4 ||p(.,t)||o,2^ t + 

+ |X n;t \^,t| 1/4 ||p„(.,t)|| 0> 2,K tl , t ) + C||VT n ||g i2i0 , 

(5.130) 

where C = C(0, Z\, . . . , Zj, e, «i, k®). Then, integrating (|5.130[) in (0,T), we 
obtain 

/ / |V/„| 2 dxdr< 
Jo Jo 

< C(B, + ||*||i A eo)(B, + ||po||o,2,jf ) / (I^W| 1/4 + |#»,AA,,t| 1/4 )dt+ 

Jo 

J T J / T \ 

+E/ ii(M*-r„(M:j)(.,i)iio,2,o^+^L2 /2 / i^n,A^r /2 ^ + 

4=1 J 4=1 \ J0 / 

J / T \ ^ 2 T 

+ Y, L 2 2 [l iKMn^dt) +^ ||(p-Pn)(.,i)||0,2,A n , t ^+ 

+ CT||VT n ||g )ai0 , 

(5.131) 

using Schwarz's inequality, ((2~5j) , (|4.112[) and dUSDJ). Now, denoting by Q(t), Q„(t) 
the afhnc isometries such that K t = Q(t)Ko, K n; t = Q n (t)Ko, we have 

||On(.)-0(-)llo ) oo,(o,T)->0, ra->oo (5.132) 
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using (|5.124l) (see [H]). As a consequence, 



\\{p- P n)(;t)\\l2.A nt < / MQn(*)e,*)-Pn(Qn(*)£,i))X 

JK Q 

(po(Q- l (t)Q n (t)t,t) - poAOfdt 



< 2 / (poiQ- 1 (t)Q n (t)0 - PoiQ-WQm) 2 ^- 

JKq 



(5.133) 



+ 2 / (p (0 - po,»(0) -> 0, n -> oo, 

using L 2 -continuity and the strong convergence ||po — Po,n||o,2,_fs" — 0, n — s> oo. 

Moreover, from (|5.132p . we have \K n j\A n ^ \ — > 0, \K t \A n ^\ — > as n — > oo. 
Then (15325)) follows from (j5T3Tj) . (j5T33l . (j5T24| . (j5T27l) and the dominated 
convergence theorem. □ 



Lemma 5.4. For each i 6 {1,...,J}, the junction M* , satisfies H3.38\) , con- 
sidering u and ^ given in \5.12J$ and Lemma \5.3\ respectively. 

Proof. For each c 6 iJ 1 (fi), we consider its extension <f€ iJ 1 (C x (0,T)) and a 
sequence G C^"* 1 ^ x [Q, T ]) such that 

-?lli,2,Ox(o,T) -> 0, n -s- oo. (5.134) 

Let us set ? n := ? n \^ T , where fl n ,T = [J {t} X f2 ni i and, as before, f2 ni j 

te(o,T) 

denote the fluid domains related to lZ n (u n ). We want to pass to the limit in 

Kn^ndx I (t) + 



f f <; n (n n (u n )-VAf* n )dxdT- 
Jo Jn n , T 

di Jo L ( VA/ ^' n + i^ K * nVm ") ■ v?n dxdr+ 

J\f* n d t <; n dxdr = / M,o,n <To,nrfx I (£) 
\Jn / 



JQ, 



for a. e. t € (0,T). 

First we observe that 



■Ni ,0,n Co.n^ x ' 



Jn 



0, n^-oo, (5.135) 



follows directly from (|5.134j) and from the strong convergence M,o,n — > A/i.o m 

Now, from (|5.124p . up to a extraction of a subsequence, we have 
Af* n (.,t)^JS*(.,t) in L 2 (0), 
a. e. i € (0,T), Vi G {1, . . . , J}. Then, using the notation of the Lemma [5.31 
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and Holder's inequality, we have 



f A/^dx - [ MfsdK < [ \K*n ~ W | |?„ |dx + f |fe - ?| |A? |dx+ 

|A£„c„|<ix + / pV?c|dx 
< - JV?(.,t)||o, 2 ,o||^(-,*)llo,2,o + IW(-,*)||o,2,o||(&» - ?)(.,t)||o,a,o+ 

+ |^A^,t| 1/6 |l^:„(-,«)||0,2,K t ||?n(-,i)||0,3,O+ 

+ |^'n,A^n,i| 1/6 ||^?(.,t)||0,2,K„, t ||?'(.,t)||0,3,O "> «, U -> OO, 

(5.136) 



for a. e. t G (0,T). Indeed, from (|5.134p . ||?„||i,2,Ox(o,r) < ||?l|i,2,Ox(o,r) for n 
sufficiently large. As a consequence, using (|5.124l) and (|5.134l) the first two terms 
above tends to zero as n — > oo. As in Lemma 15 .31 |ATA i A n)t |, \K n ^\A n ^\ — > 0, 
n oo; then from Sobolev embedding and (I4.112p . we obtain the convergence 
of the third and fourth terms above. 

From (|5.113|) , we have 



VAf* n ^VAT* weakly in L 2 ((0,T) x Of 



(5.137) 



Vi e {I,..., J}. For a.e.t£ (0,T), 



[* [ Sn (ftn(u„) • VA/^Jdxdr - / / ? (u ■ VjV*)dxdr 
JO J£l„, x JO Jf2 T 

< / / |VA^: n ||f„-?l|^ n (u n )|dxrfr+ 
Jo Jsi„, T 

/ / ?u ■ VA/f dxdr - f f ?u-VACdxdr 

Jo Ja„ iX Jo Jn T 

f f ?ll n (u n )-VN* n dxdT- f f ?u- VA/^dxdr 

JO JQ„,t JO JQn.T 



— > 0, n — >■ oo. 
(5.138) 



In fact, Holder and Poincare inequality, Sobolev embedding, (|4.109p . (|4. 1 12[) . 
(|5.134p and (|4.64j) furnish us with the estimate 



/ / |V/v7J?„-?l|ft n (u„)|dxd r 
Jo Jn n t 



< / l|VA^* n || ,2,o||?n-?l|o,3,o||^n(u„)||o,6,orfr 
JO 



<C(0 > T)||g l -?l| 1 , a>0x( o,T)| / ||VA^ n ||g ia>0 dr 



1/2 



1/2 



||V^„(u n )||2 



2 : 



< C(0,T,i 1 ,i 2 ,C)[[? n -?1| 1)2 ,c x( o,T)^0, n^oo. 
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Also, we have, 



/ f ?u • VJV* n dxdT - [ [ qu-VAf*dxdr 
Jo Jn n<r Jo Jn T 

n x (o,t) (2h„ T ~ ln T )?u • VAf* n dxdr 
j 

[ T [ l (0tt) l n ^u-(VAr* n -VAf*)dxdT 
Jo Jo 



The second term in the right-hand side above tends to zero as n — > oo, using 
([5137]) and the fact that l (0}t) X UT Zu E L 2 (0 x (0,T)) 3 . From (^T32j) . l Un T - 
1-Q, t — > a. e. in O x (0,T), with n — > oo. Then, 



nI(o,t) (Ia„ T - Za T )?u • X7J\f* n d-xdT 
j 



< / ll^.x - I oJ|o,8,o||VA/; H ; il ||o,2,o||?l|o,16/3,ol|u||o,16/3,C'^ 

< C||?lli,2,e)x(o,T) / l|2«„, T -2n T ||o,8,o||VAi* n ||o,2,o||u||o,i6/3,o^''" 

\ 1/2 







< C||?lll,2,Ox(0,T) / ll u llo46/3,C>ll :Z: f2„, T -^110,8,0^ 











\\VK*Jh,odT 



<cm\ 



l,2,Ox(0,T) 



1/2 



10,16/3 



,o\\ x ^,r ~ In T \\o,s.,o dT -> °> n oo, 



where C = C(0,T) and we have used the dominated convergence theorem. 
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Also, from Holder's inequality, Sobolev embedding, (|4.112|) . (I4.109j) and (I5.124p . 

n?K n (u n ) ■ VAf* n dxdr - [ [ ?u ■ VAf* n dxdT 
-i n ,T Jo Ju„, T 

< I I mn n (u n ) - u\\VM* n \d X dr 
Jo Jo 



< 



fl^K) - ulV^lv^J 1 / 2 ^^^ II \K n (u n ) - ul^lv^l^Ho^odr 



f T 

< / ||?l|o,6,o||^n(u„) ~ U||y« o ||V^ n || 0l 2,o||^(u„) - U^ 2 Q dT 
JO 

< C(O,T)||?l| 1)2iOx{0)T) ^jf \\K n (u n ) - u|| 2 / 2 3 ||VA^:J|^ 2 3 dr^ 

X ^ ll^nK) -u|| 2 2 dT^J 

<C(O,T)\\?\\ h2!Oxi0!T) (^ ||7^(u n ) - uH^odr") x 



x (j* \\K n (u n ) - uWl ^dr^ (j* IIVA/^II 2 , 



1/2 



2,0 



dT 



<C(0,T,L 2 ,L 1 ,C)(j^ ||^„(u„)-u|| 2 2:0 dr N J -> 0, n ->• 00. 



The convergences 



/ / VA/"* • V<?„ dxdr - [ [ VA"* • V<? dxdr 
'0 Jn„ iT Jo Jsi T 



->0, 



A7„ d t <,n dx- I I N* d t c, dxdr 



Jn„ 



Jn T 



(5.139) 



— > 0, n — > oc 



for a. e. t £ (0,T), follows from analogous arguments. 
Now, for a. e. t 6 (0,T), 



/ / • V<;„ dxdr - f [ N*V^ • rfxdr 

Jo Jq„,t Jo Jo T 

< [ T [ AA*„|V[^ n ]"||V^-V^xdr+ / / ACJV[^„] n - W||V^|dxdr- 

JO JS„ iT JO JB„ |T 

+ / T / |A^-A7||V^||V^xdT+ 
Jo Jb„. t 

|A^„||V^„] n ||V^| dxdT+ 



10 JK T \A n , T 
pT 



|A/?||VV>||V?| dxdr -> 0, 



Jif n , x \yl„ 



(5.140) 
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In fact, from Holder's inequality, Sobolev embedding, P~B2|) , (|4.110|) . (|4.112|) 
and (|5.134l) we have 



/ / |^ n ||V[Vn] B ||Vfe-Vf|&dr 

JO JB„, T 

< / ||V£ - V?l|o,2,o||A^*„|| ,4,o||V[^]"||o,4,n„ T dr 
Jo 

<C(O,T)||^-?1| li2)Ox(0iT) x 

/ T \ ^ 2 / T \ 

X / H^nlllAO*] / l|V[Vn] n ||2,4, nBT drJ 



< C(0,B,,L 2 ,C,T)|K~ -?1|i, 2i0x (o,t) ->0, »->oo 



and 



/•T 



■A/^|V^„] n -V^||V<r|dxdr 

y |V51 2 dxj dr 



T / \ / \ I/ 2 

-I (X iM:,j 2 iv[^r-v^i 2 dxj 

rp / \ 3/8 

<C(0,T)||?1| 1 , aiOx( o, T) / [/ |^ n | 8 / 3 |V[^] n -V^| 2 dx x 



\JB„ 

1/8 



x (J WbPnF - V^| 2 dx^ dr 
x f ll^:J|o,i6/3,ol|V[^]" - V^|lof 4 %„ ,J|VhM n - ^Wl'lodr 



T \ 1/2 



<C(p,T)m\w^ T) ^ ||A^:J| 2 , 2i0 drJ x 

\ 3/8 / ,T 

x(/ HVW-VVHgoodr / ||Vty„] n -V^||g AO dT 



o 

T v 3/8 / T v 1/8 



II 



<C(0,T,B.,L2)\\s\\i,2,oxQ,r)U l|V[^„] n - V^l^dr j -> 0, n -> oo 

(5.141) 

Analogously, 

/ / -A/J*||V^||V?|<ixdr -> 0, n oo. 

Jo Js„. T 
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Finally, we have, 

f \Af*\\ViP\\V<;\dxdT 

JK n , r \A n , T 

< [ |^„,t\A„ >t | 1 / 12 ||A^*||o,6,o||VV'||o J 4,K„, T ||V?l|o,2,odT 



1/2 / „ T \ 1/2 

\\K*\\h,odT) 



<C(0 ) T,B»)||?1| 1 , 2 , 0x( o,t) (j a \K nit \An,t\ 1/6 dt^j ^ 
<C(O,T,B*,L 2 ,CW\\i.2,ox(0,T)( [ \K n , t \An,t\ 1/6 dt) -> 0, n -»• oo, 



(5.142) 

using the dominated convergence theorem. Also, 

/ |-A/£J|VfoM n ||Vc n |dxdr 

JK T \A n , T 

< f [ l^nllVW - V^||V^| dxdr+ 

JO JK T \A n , T 

+ [ I W* n W\\V<; n \ dxdr^Q, n^oo, 

JO JK T \A n<T 

using similar arguments as in (|5.142l) and (|5.14ip . The result follows from 

(ESS), HSIT3S}, (jSTTSgi) and flOU). □ 

Lemma 5.5. The function u satisfies {3.3$ , considering ip, A/i and /i given in 
15.3]) . {5.12$ and {5.11$ , respectively. 



Proof. The proof follows from (|5.124p and from the arguments in [15] , observing 
that: given w e S and w„ € S n such that L T ||w n — w||f 2 d^ t with n — > oo, 
we obtain 

F„ • w n dr — / F • wdr — > 0, ri — S> oc 
./o 

for F„ = -eEtA/^Vi, F = -e£f =1 Z^K*^, this follows 
from analogous estimates as in Lemma 15.41 □ 

We have obtained T > and a weak solution for the problem in (0, T). From 
the bound 1)4. 109)1 . Q(t) is continuous in [0,T] (see [14]). so that 

lim dist( K t ,dO) > 7 > 

and we can iterate the existence result in order to obtain Theorem 13. II 
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